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CHAPTER 1 
GENERAL 



1 -1 . Purpose and scope. 
This manual contains criteria and 
methods for calculating the depths of 
freeze and thaw in soils, with considera- 
tion of the effects of other adjacent 
materials, for the design of military 
facilities in seasonal frost, arctic and 
subarctic regions. The contents are 
applicable to both Army and Air Force 
construction in arctic, subarctic and 
seasonal frost areas. The data pre- 
sented in this manual relate to arctic 
and subarctic facility design presented 
in the other manuals of the Arctic and 
Subarctic Construction series. 

1 -2. References and symbols. 
Appendix C lists the references for 
this manual; appendix A contains a list 
of symbols. 

1 -3. Background. 

a. The depths to which soils may 
freeze and thaw is very important in 
the design of pavements, structures 
and utilities in areas of seasonal frost 
and permafrost. Methods of calculating 
such depths, based on heat-transfer 
principles, are presented here. For 
derivation of basic equations, and the 



underlying theory, see appendix C and 
the bibliography. 

b. Heat transfer in soils involving 
phase change of pore water is an ex- 
tremely complex process and many 
problems defy rigorous mathematical 
treatment. The methods presented here 
are simplified procedures developed 
for the solution of engineering design 
problems. 

c. Several assumptions have been 
made in developing practical methods 
of calculating depths of freeze or thaw 
in soils. It is assumed that each layer of 
material is homogeneous and isotropic, 
and that the average thermal proper- 
ties of frozen and unfrozen soils are 
applicable. Unless specific data are 
available, it is also assumed that all soil 
water is converted to ice, or all ice is 
converted to water, at a temperature of 
32°F. This latter assumption is sub- 
stantially correct for coarse-grained 
soils but only partially true for fine- 
grained soils. 

d. The services of the U.S. Army Cold 
Regions Research and Engineering 
Laboratory (US ACRREL), Hanover, New 
Hampshire, are available to assist in 
the development of solutions for heat- 
flow problems in soils. 
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CHAPTER 2 
DEFINITIONS AND THERMAL PROPERTIES 



2-1. Definitions. 

Definitions of certain specialized terms 
applicab le to aj^cjtic and subarctic 
r egio ns are contai ned in TM 5- 8 52/ 
AFR "8849, Volume 1. Following are 
additional terms used specifically in 
heat-transfer calculations. 

a. Thermal conductivity, K. The 
quantity of heat flow in a unit time 
through a unit area of a substance 
caused by a unit thermal gradient. 

b. Specific heat, c. The quantity of 
heat absorbed (or given up) by a unit 
weight of a substance when its tempera- 
ture is increased (or decreased) by 1 
degree Fahrenheit (°F) divided by the 
quantity of heat absorbed (or given up) 
by a unit weight of water when its 
temperature is increased (or decreased) 
by 1°F. 

c. Volumetric heat capacity, C. The 
quantity of heat required to change 
the temperature of a unit volume by 
1°F. 

— For unfrozen soils, 



yd (c + 1.0 



lOO 



— For frozen soils, 



C f = yd (c + 0.5 



lOO 



(eq2-l) 



(eq2-2) 



-Average values for most soils, 



C avg " Yd <° + °" 75 



lOO 



(eq 2-3) 



where c = specific heat of the soil solids 
(0.17 for most soils) 

y A = dry unit weight of soil 
w = water content of soil in per- 
cent of dry weight. 
d. Volumetric latent heat of fusion, 
L. The quantity of heat required to 
melt the ice (or freeze the water) in a 
unit volume of soil without a change in 
temperature— in British thermal units 
(Btu) per cubic foot (ft 3 ): 



L = 144 



Xd 



lOO 



(eq2-4) 



e. thermal resistance, R. The recipro- 
cal of the time rate of heat flow through 
a unit area of a soil layer of given thick- 
ness d per unit temperature difference: 



Ri 



(eq 2-5) 



K 



/. Thermal diffusivity, a. An indica- 
tor of ho w easily a material will undergo 
temperature change: 



K 



g. Thermal ratio, a . 



<eq2-6) 



(eq2-7) 



where 

v = absolute value of the difference 
between the mean annual tem- 
perature below the ground 
surface and 32°F. 

v a = one of two possible meanings, 
depending on the problem 
being studied: 

(1) v s = nF/t (or nl/t) 
where 

n = conversion factor from air 
index to surface index 

F= air freezing index 

I = air thawing index 

t = length of freezing (or thawing) 
season. 

(2) v s = absolute value of the dif- 

ference between the mean 
annual ground surface temp- 
erature and 32°F. 

(3) In the first case, v s is useful for 
computing the seasonal depth 
of freeze or thaw. In the second 
case, it is useful in computing 
multiyear freeze or thaw 
depths that may develop as a 
result of some long-term 
change in the heat balance at 
the ground surface. 

h. Fusion parameter, //• 

2-1 
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u= v g (eq2-8) 

L, 

where v s has the two possible meanings 
noted above. 

t "Lambda" coefficient, A- A factor 
allowing for heat capacity and initial 
temperature of the ground (see fig. 3-1). 



a = f(a,/Lt).* 



(eq 2-9) 



j. Thermal regime. The temperature 
pattern existing in a soil body in rela- 
tion to seasonal variations. 

k. British thermal unit, Btu. The 
quantity of heat required to raise the 
temperature ofl pound (lb) of water 1°P 
at about 40°F. 



2-2. Thermal properties of soils and other 
construction materials. 

a. The basic thermal properties of 
soils and other construction materials 
used to calculate depths of freeze and 
thaw are specific heat, thermal con- 
ductivity and volumetric latent heat of 
fusion. Other terms used in heat-flow 
calculations are derived from these 
data and the elements of weight, length, 
temperature and time. 

b. the specific heat of most dry soils 
near the freezing point may be assumed 
to be constant at the value of 0.17 Btu/lb 
°P. Specific heats of various materials 
are given in table 2-1; see the ASHRAE 
Guide and Data Booh of the American 



Table 2-1. Specific heat values of various materials 
(U.S. Army Corps of Engineers}. 



Material 



Aluminum 
Asbestos fibers 
Concrete (avg. stone) 
Concrete (dams) 
Copper 
Corkboard 

Cork, granulated 
Fiberglas board 

Foamglas 

Glass block, expanded 

Glass sheets 

Glass wool 

Ice 

Iron (alpha) 

Masonry 

Mineral wool 

Perlite, expanded 

Polystyrene, cellular foam 

Polyurethane foam 

Sawdust 

Snow 

Steel 

Straw 

Water 

Woods (avg.) 

Woods fiberboard 



, Temperature 


Specific heat 


TF) 


(Btu/lb °F) 


-27.4 


0.20 


— 


0.25 


— 


0.20 


~ 


0.22 


44.6 


0.20 


91 


0.43 


-19 


0.29 


— 


0.42 


111 


0.24 


-22 


0.19 


-20 


0.16 


112 


0.18 


— 


0.20 


— 


0.16 


32 


0.48 


44.6 


0.11 


— 


0.22 


— 


0.22 


— 


0.22 


— 


0.27 


— 


0.25 


— 


0.60 


— 


0.50 


_. 


0.12 


— 


0.35 


— 


1.00 


68 


0.33 


148 


0.34 



Specific heat values shown to nearest 0.01. Average values listed where temperature is not shown. 
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Society of Heating and Air Condition- 
ing Engineers for the specific heat 
values of common materials. 

c. The thermal conductivity of soils 
is dependent upon a number of factors: 
density; moisture content; particle 
shape; temperature; solid, liquid and 
vapor constituents; and the state of the 
pore water, whether frozen or unfrozen. 
Average values, expressed in Btu/ft 
hour °F, for frozen and unfrozen granu- 
lar soils, silts and clays should be read 
from figures 2-1 through 2-4. The 
charts for sands and gravels are applic- 
able when the silt and clay content 
together make up less than 20% of the 
soil solids. The charts for silt and clay 
are applicable when that fraction is at 
least 50%. For intermediate silt-clay 
fractions, it is recommended that the 
simple average of the values for the 
two sets of charts be used. In all cases, 
the error in the thermal conductivity 
estimates may be ±25%, and even 
higher when the percentage of quartz 
grains in the soil is exceptionally high 
or low. Figures 2-5and2-6 presentesti- 
mates of the average thermal conduc- 
tivity of frozen and unfrozen peat. An 
excellent source of data for dry con- 
struction materials is the ASHRAE 
Guide and Data Book. Thermal con- 
ductivity values for a number of com- 
mon construction materials are listed 
in table 2-2. 

d. The latent heat of fusion is the 
amount of heat required to cause a 
phase change in soil moisture. This 
amount of heat does not change the 
temperature of the system when 
freezing or thawing takes place. The 
gravimetric latent heat of fusion of 
water is assumed to be 144 Btu/lb. The 
amountofheat energy required to con- 
vert 1 ft of water to ice is (144 X 62.4 =) 
9000 Btu/ft 3 and to change 1 ft 3 of ice to 
water is (144 X 0.917 X 62.4 =) 8240 
Btu/ft 3 . (Note: The density of water 62.4 
lb/ft 3 and that of pure ice is 57.2 lb.ft 3 ). 

e. Figures 2-7 and 2-8 may be used 
to determine the average volumetric 
heat capacity and volumetric latent 
heat of fusion, respectively, of moist 
soils. 

/. The following ea;ampZe illustrates 
the significance of latent heat of fusion 
relative to th,e volumetric heat capacity 
for a moist soil. Assume a soil having a 
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dry unit weightof 120 lb/ft 3 and a water 
content of 15 percent. Its volumetric 
latent heat of fusion, L, is (144 x 120 X 
0.15 =\ 2592 Btu/ft 3 and its average volu - 
metric heat capacity, C, is (120 [0.17 + 
0.75 X 0.15] =) 33.9 Btu/ft 3 °F. The quan- 
tity of heat required to change the 
phase of pore water in 1 ft 3 of this soil at 
32°F is the same as that required to 
cause a temperature change of 
(2592/33.9 =) 76.4°F when a phase 
change is not involved. 

2-3. Fundamental considerations. 

a. Theoretical basis. The freezing 
or thawing of soils is the result of re- 
moving or adding heat to an existing 
soil mass. The movement of heat is 
always in the direction of lower temper- 
ature. The time rate of change of heat 
content depends on the temperature 
differential in the direction of heat flow 
and on the thermal properties of the 
soil. 

b. Physical factors and data re- 
quired. Calculation of the depth of 
freeze or thaw is based on knowledge 
of the physical and thermal properties 
of the soil in the profile, the existing 
thermal regime, and the nature and 
duration of boundary conditions 
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Figure 2-1. Average thermal conductivity tor sands and 
gravels, frozen. 
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(U.S. Army Corps of Engineers) 

Figure 2-2. A verage thermal conductivity for sands and gravels, unfrozen. 



190 



Y d 




90 


120 


100X 


saturated 


50 




N/A 


N/A 




1.20 


40 




1.00 


N/A 




1.23 


30 






0.98 




1.26 


20 
pure 


Ice 




0.61 




1.26 
1.26 



Thermal conductivity K is expressed In Btu per hour per square foot per 
unit thermal gradient in °F per foot: 

Dashed line represents extrapolation. 

(U.S. Army Corps of Engineers) 

Figure 2-3. Average thermal conductivity for silt and clay 
soils, frozen. 




10 20 SO 

MOISTURE CONTENT, Ptrctnt 
(U.S. Army Corps of Engineers) 

Figure 2-4. Average thermal conductivity for silt and clay soils, unfrozen. 
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Figure 2-5. Average thermal conductivity for peat, frozen. 
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Figure 2-6. Average thermal conductivity for peat, unfrozen. 
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causing a change in the thermal regime. 
Data pertinent to the soil profile in- 
clude grain-size distribution, classifica- 
tion, density, water or ice content, and 
temperature of each soil stratum. 
Knowledge of the thermal properties 
of all the materials in the heat flow 
path is also required. Measured or as- 
sumed temperatures within the soil 
mass determine the initial conditions. 
If surface temperatures can be assumed 
to be spatially uniform and the thermal 
influence of any buried structures can 
be considered negligible, one-dimen- 



sional heat flow may be assumed, there- 
by simplifying the problem and its 
solution. 

2-4. Freezing and thawing indexes. 

a. Physical concept and quantita- 
tive measurement. The penetration of 
freezing or thawing temperatures into 
soil partly depends on the magnitude 
and duration of the temperature dif- 
ferential at the air-ground interface. 
The magnitude of the temperature dif- 
ferential is expressed as the number of 
degrees that the temperature in the air 

2-5 
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Table 2-2. Thermal properties of construction materials 
(U.S. Army Corps of Engineers). 



Type of 




material 


Description 


Asphalt 


Mix with 6% by 


paving „ 


weight cut- 


mixture 


back asphalt 


Concrete 


With sand and 




gravel or stone 




aggregate (oven- 




dried) 




With sand and 




gravel or stone 




aggregate (not 




dried) 




With lightweight 




aggregates, in- 




cluding expanded 




shale, clay or 




slate; expanded 




slags; cinders; 




pumice; perlite; 




vermiculite; also 




cellular con- 




cretes. 


Wood 


Maple, oak and 




similar hard- 




woods 




Fir, pine and 




similar soft- 




woods 


Building 


Asbestos- 


boards 


cement board 




Plywood 




Wood fiberboard, 




laminated or 




homogeneous 




Wood fiber- 




hardboard type 


Blanket 


Mineral wool, 


and batt 


fibrous form, 


insula- 


processed from 


tion 


rock, slag, or 




glass 




Wood fiber 


Board 


Cellular glass 


and 


Corkboard (with- 


slab 


out added binder) 


insula- 


Glass fiber 


tion 


Wood or cane fiber- 




interior finish 




(plank, tile, lath) 




Expanded polystyrene 




Expanded ureaformal- 




dehyde 




Expanded perlite 




Polyurethane foam 




Mineral wool with 




resin binder 




Mineral wool with 


2-6 


asphalt binder 



Unit weight 
(lb/ft 3 ) 



138 



140 



140 
120 

100 
BO 
60 
40 
30 
20 



45 



32 



1.5-4.0 



conductivity 
(Btu/ft 2 -br-°F 
per in.) 



10.3 



9.0 



12.0 
5.2 

3.6 

2.5 

1.7 

1.15 

0.90 

0.70 



1.10 



0.80 



0.27 



3.2-3.6 


0.25 


9.0 


t0.39 


6.5-8.0 


t0.27 


9.5-11.0 


0.25 


15.0 


0.35 


1.6 


0.29 


1.0 


0.25 


9.5-11.5 


0.34 


1.5-3.0 


0.17 


15.0 


**0.28 


15.0 


**0.31 



K 
conductivity 
(Btultt- 
hr-°F) 



0.86 



0.75 



1.00 
0.43 

0.30 

0.21 

0.14 

0.096 

0.075 

0.058 



0.092 



0.067 



120 


4.0 


0.33 


34 


0.80 


0.067 


26,33 


0.42,0.55 


0.035,0.046 


65 


1.40 


0.12 



0.022 

0.021 
0.032 

0.022 
0.021 

0.029 

0.024 



0.023 
0.026 



*TM 5-852-6/AFR 88-19, Volume 6 



Table 2-2. Thermal properties of construction materials, Continued 









conductivity 


conductivity 




Type of 




Unit weight 


(Btu/ft z -br-°F 


(Btu/ft- 




material 


Description 
Cork, granulated 


(lb/ft 3 ) 
5-12 


per in.) 
0.25-0.36 


hr-°F) 




Loose 







fill 


Expanded perlite 


3-4 


0.28 


— 




insula- 


Mineral wool 










tion 


(glass, slag, or rock) 


2-5 


0.30 


0.025 






Sawdust or shavings 


8-15 


0.45 


0.037 






Straw 


7-8 


0.32 








Vermiculite (ex- 












panded) 


7.0-8.2 


0.48 


0.040 






Wood fiber: redwood, 












hemlock, or fir 


2.0-3.5 


0.30 


0.025 




Miscel- 


Aluminum 


168 


1416 


118 




laneous 


Copper 


549 


2640 


220 






Ductile iron 


468 


360 


30 






Glass 


164 


5.5 


0.46 






Ice 


57 


"15.4 


1.28 






Snow, new, loose 


5.3 


0.6 


0.05 






Snow on ground 


18.7 


1.56 


0.13 






Snow, drifted and 












compacted 


31.2 


4.8 


0.40 






Steel 


487 


310 


25.8 






Water, average 


62.4 


4.2 


0.35 






Values fork are for dry building materials at a mean temperature of 75° F except as noted; wet conditions will adversely affect values of 
many of these materials. 
tMean temperature of 60°F. 
"Mean temperature of 32° F. 



or at the ground surface is above (posi- 
tive) or below (negative) 32°F, the as- 
sumed freezing point of water. The 
duration is expressed in days. 

b. Air freezing index (F) and air 
thawing index (I). The air freezing and 
air thawing indexes, as defined in TM 
5-852-1/ AFR 88-19, Volume 1, may be 
determined by the following methods. 
(1) Summation of degree-days of 
freeze and thaw from average daily 
temperatures. If T x is the maximum 
daily air temperature and T 2 is the 
minimum daily air temperature, the 
average daily air temperature T, may 
be taken as 1/2 (T 1 + T z ), and the number 
of degree-days for the day is (T - 32) . 
The summation of the degree-days for 
a freezing or thawing season gives the 
air freezing or air thawing index. 
Table 2-3 illustrates the method used 
to obtain the summation of degree-day s 
for a 1- week period, assuming that -456 
degree-days had been accumulated 
since freezing began. An average dally 
temperature based on hourly tempera- 
tures would be slightly more accurate, 
but such precision is not usually war- 
ranted. The negative sign, indicating 



freezing degree-days, is usually 
omitted. 

(2) Calculation fromaverage month- 
ly temperatures. The freezing or 
thawing index may be calculated by 
determining the area between the 32°F 
line and the curve of average monthly 
temperature and time, taken over the 
appropriate season. The area may be 
determined by planimeter or a simple 
approximation rule (Simpson's rule, 
midordinate rule, etc.). The areas are 
expressed in units of degree-days, re- 
sulting in a summation of degree-days 
or a freezing or thawing index. For an 
example refer to figure 2-9, a plot of 
the monthly average temperatures at 
Fairbanks, Alaska, from September 
1949 to October 1950. Determination of 
areas by planimeter gave a freezing 
index of -5240 degree-days and a 
thawing index of +3420 degree-days. 
The use of Simpson's rule gave a freez- 
ing index of -5390 and a thawing index 
of +3460 degree -days. Either pair of 
indices is adequate for computations. 

c. Surface-freezing and surface- 
thawing indexes. For determining the 
heat flow within the soil, it is neces- 
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Figure 2-7. A verage volumetric heat capacity for soils. 



sary to determine or estimate the temp- 
erature condition at the ground sur- 
face. Since air temperatures are gen- 
erally available and surface tempera- 
tures are not, a correlation between 
these temperatures helps establish the 
thermal boundary condition at the 
ground surface. The combined effects 
of radiative, convective and conductive 
heat exchange at the air-ground inter- 
face often must be considered in deter- 
mining surface temperature. 

d. Correlation of air and surface 
indexes. No simple correlation exists 
between air and surface indexes. The 
difference between air and surface 
temperatures at any specific time is 
influenced by latitude, cloud cover, 
time of year, time of day, atmospheric 
humidity and stability, wind speed, 
snow cover and ground surface char- 

2-8 



acteristics, and subsurface thermal 
properties. Heat balance algorithms 
that approximate many of these inter- 
relationships exist but they are often 
unwieldly to use and the inputs are 
often difficult to characterize. It is 
recommended that the ratio of surface 
index to air index, designated as the 
"n-factor," be used to represent a 
monthly or seasonal correlation. Re- 
liable determination of the n-factor for 
a specific location requires concurrent 
observations of air and surface tem- 
peratures throughout a number of com- 
plete freezing and thawing seasons 
plus anticipation of future changes to 
conditions existing during the period 
of measurement. Such determination 
is often not feasible, so n-factors must 
generally be estimated conservatively 
from n-factors tabulated for other, pref- 
erably similar, sites. 
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Figure 2-8. Volumetric latent heat for soils. 



Table 2-3. Calculation of cumulative degree-days 
(U.S. Army Corps of Engineers). 





Temperature (° F) 






Degree-days 


Cumulative 


Day 


Maximum 


Minimum 


Average 
15 


per day 
-17 


degree-days 


1 


29 


1 


-473* 


2 


9 


-11 


-1 


-33 


-506 


3 


10 


-8 


1 


-31 


-537 


4 


15 


-1 


7 


-25 


-562 


5 


30 


16 


23 


-9 


-571 


6 


38 


30 


34 


2 


-569 


7 


30 


18 


24 


-8 


-577 



"Prior accumulation of -456 degree-days assumed. 
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Figure 2-9. Average monthly temperatures versus time at Fairbanks. Alaska. 



(1) Freezing conditions. The n- 
factor is very significant in analytical 
ground studies. It generally increases 
with wind speed. Snow cover reflects a 
large part of incoming solar radiation 
resulting in higher freezing indexes 
at the snow surface, but its insulating 
effect can greatly reduce the freezing 
index at the ground surface. The effects 
of turf or an organic ground cover on 
the heat flow processes at the air- 
ground interface are extremely vari- 
able and difficult to evaluate. On the 

2-10 



basis of observations and studies made 
to date, the n-factors given for average 
conditions in table 2-4 should be used 
to convert the air freezing index to the 
surface freezing index in the absence 
of specific measurements at the site of 
planned construction. 

(2) Thawing conditions. The n- 
factor for thawing conditions is af- 
fected by the same factors as those for 
freezing conditions. It is the ratio of 
surface degree-days of thaw (degrees 
above 32°F) to air degree-days of thaw. 



Incoming shortwave radiation may in- 
troduce heat into the surface to an 
extent that the surface becomes a 
source of heat conducted not only down- 
ward but upward into the air. In such a 
case the n-factor may become signifi- 
cantly larger than 1.0. The effect of 
latitude is not particularly significant 
in arctic and subarctic areas, but con- 
sideration should be given to the effect 
of wind speed. Recommended curves 
for n-factors versus wind speed for 
portland-cement-concrete and bitu- 
minous-concrete pavements are shown 
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in figure 2-10 and are based on field 
studies conducted in Alaska and 
Greenland. The n-factors given for 
average conditions in table 2-4 should 
be used to convert air thawing indexes 
to surface thawing indexes in the ab- 
sence of specific measurements at the 
planned construction sites. 

e. Design indexes. For design of per- 
manent pavements, the design freezing 
(or thawing) index should be the aver- 
age air freezing (or thawing) index of 
the three coldest winters (or warmest 
summers) in the latest 30 years of 



Table 2-4. n-factors for freeze and thaw (ratio of surface index to air index) 
(U.S. Army Corps of Engineers). 



Type of surface' 



For freezing 
conditions 



For thawing 
conditions 



Snow surface 

Portland-cement concrete 

Bituminous pavement 

Bare soil 

Shaded surface 

Turf 

Tree-covered 



1.0 

0.75 

0.7 

0.7 

0.9 

0.5 

0.3" 



1.5 

1.6to2f 

1.4to2f 

1.0 

0.8 

0.4 



Surface exposed directly to sun or air without any overlying dust, soil, snow or ice, except as noted otherwise, and with no building 
heat involved, 
t Use lowest value except in extremely high latitudes or at high elevations where a major portion of summer heating is from solar 
radiation. 
Data from Fairbanks, Alaska, for single season with snow cover permitted to accumulate naturally. 
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Figure 2-10. Relationship between wind speed and n-factor during thawing 

season. 
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record. If 30 years of record are not 
available, the air freezing (or thawing) 
index for the coldest winter (or 
warmest summer) in the latest 10-year 
period may be used. For design of 
foundations for average permanent 
structures, the design freezing (or 
thawing) index should be computed 
for the coldest (or warmest) winter in 
30 years of record or should be esti- 
mated to correspond with this fre- 
quency if the number of years of record 
is limited. Periods of record used 
should be the latest available. To avoid 
the necessity for adopting a new and 
only slightly different freezing (or 
thawing) index each year, the design 
index at a site with continuing con- 



struction need not be changed more 
often than once in 5 years unless the 
more recent temperature records in- 
dicate a significant change. The distri- 
bution of design freezing and thawing 
index values in North America is 
presented in TM 5-852-1/AFR 88-19, 
Volume 1. The relatively linear relation- 
ship between recorded maximum in- 
dexes and mean freezing and thawing 
indexes shown in figures 2-11 and 2-12 
may be used in conjunction with dis- 
tribution of mean and freezing and 
thawing indexes in TM 5-852-1/AFR 
88-19, Volume 1 to determine the design 
index values for arctic and subarctic 
regions. 
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Figure 2-11. Relationship between mean freezing index and maximum freezing index for 10 
years of record, 1953-1962 (arctic and subarctic regions). 
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Figure 2-12. Relationship between mean thawing index and maximum thawing index for 10 years ot 
record, 1953-1962 (arctic and subarctic regions). 
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CHAPTER 3 
ONE-DIMENSIONAL LINEAR AND PERIODIC HEAT FLOW 



3-1. Thermal regime. 

The seasonal depths of frost and thaw 
penetration in soils depends upon the 
thermal properties of the soil mass, the 
surface temperature (upper boundary > 
condition) and the thermal regime of 
the soil at the start of the freezing or 
thawing season. Many methods are 
available to estimate frost and thaw 
penetration depths and surface tem- 
peratures. Some of these are sum- 
marized in appendix B. This chapter 
concentrates on some techniques that 
require only relatively simple hand 
calculations. For the computational 
methods discussed below, the initial 
ground temperature is assumed to uni- 
formly equal the mean annual air tem- 
perature of the particular site under 
consideration. The upper boundary 
condition is represented by the sur- 
face freezing (or thawing) index. 

3-2. Modified Berggren equation. 

a. The depth to which 32°F tem- 
peratures will penetrate into the soil 
mass is based upon the "modified" 
Berggren equation, expressed as: 



(eq3-l) 



where 

X= depth of freeze or thaw (ft) 
K= thermal conductivity of soil 

(Btu/f t hr °F) 
L= volumetric latent heat of fusion 

(Btu/ft3) 
n= conversion factor from air 

index to surface index (dimen- 

sionless) 
F= air freezing index (°F-days) 
I = air thawing index (°F-days) 





48KnF 


X - X 






L 


or 






48 K nl 



^ = coefficient that considers the 
effect of temperature changes 
in the soil mass (dimensionless). 
The x coefficient is a function of the 
freezing (or thawing) index, the mean 
annual temperature of the site, and the 
thermal properties of the soil. Freeze 
and thaw of low-moisture-con tent soils 
in the lower latitudes is greatly in- 
fluenced by this coefficient. It is deter- 
mined by two factors: the thermal ratio 
a and the fusion parameter //. These 
have been defined in paragraph 2-1. 
Figure 3-1 shows A as a function of a 
and/u. 

b. A complete development of this 
equation and a discussion of the neces- 
sary assumptions and simplifications 
made during its development are not 
presented here. A few of the more im- 
portant assumptions and some of the 
equation limitations are discussed 
below. The assumptions and limita- 
tions apply regardless of whether the 
equation is used to determine the depth 
of freeze or the depth of thaw. 

(1) Assumptions. The mathematical 
model assumes one-dimensional heat 
flow with the entire soil mass at its 
mean annual temperature (MAT) prior 
to the start of the freezing season. It 
assumes that when the freezing season 
starts, the surface temperature 
changes suddenly (as a step function) 
from the mean annual temperature to 
a temperature v s degrees below freez- 
ing and that it remains at this new 
temperature throughout the entire 
freezing season. Latent heat affects the 
model by acting as a heat sink at the 
moving frost line, and the model as- 
sumes that the soil f reezes at a tempera- 
ture of 32°F. 

(2) Limitations. The modified 
Berggren equation is able to determine 
frost penetration in areas where the 
ground below a depth of several feet 
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Figure 3-1. A coefficient in the modified Berggren formula. 
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remains permanently thawed, or to 
determine thaw penetration in areas 
where the ground below a depth of 
several feet remains permanently 
frozen. These two conditions are simi- 
lar in that the temperature gradients 
are of the same shape, although re- 
versed with respect to the 32°F line. No 
simple analytical method exists to deter- 
mine the depth of thaw in seasonal 
frost areas or the depth of freeze in 
permafrost areas, and such problems 
should be referred to HQDA (DAEN- 
ECE-G) or HQ AFESC. Numerical tech- 
niques and computer programs are 
available to solve more complex prob- 
lems. Appendix B discusses some ther- 
mal computer models for computing 
freeze and thaw depths. The modified 
Berggren equation cannot be used suc- 
cessfully to calculate penetration over 
parts of the season. The modified 
Berggren equation does not account 
for any moisture movement that may 
occur within the soil. This limitation 
would tend to result in overestimated 
frost penetration (if frost heave is sig- 
nificant) or underestimated thaw 
penetration. 

(3) Applicability. The modified 
Berggren equation is most often applic- 
able in either of two ways: to calculate 
the multi-year depth of thaw in perma- 
frost areas or to calculate the depth of 
seasonal frost penetration in seasonal 
frost areas. It is also sometimes used to 
calculate seasonal thaw penetration 
(active layer thickness) in permafrost 
areas. 

3-3. Homogeneous soils. 
The depth of freeze or thaw in one 
layer of homogeneous soil may be deter - 
mined by means of the modified 
Berggren equation. A thin bituminous 
concrete pavement will not affect the 
homogeneity of this layer in calcula- 
tions, but a portland-cement-concrete 
pavement greater than 6 inches thick 
should be treated as a multilayered 
system. In this example for homoge- 
neous soils, determine the depth of 
frost penetration into a homogeneous 
sandy silt for the following conditions: 
— Mean annua 1 temperature 
(MAT) = 37.2°F. 



—Surface freezing index (nF) = 

2500 degree-days. 
—Length of freezing season (t) 

= 160 days. 

—Soil properties: y d = 10 ° lb/ft 3 , w 

= 15%. 

The soil thermal properties are as 
follows: 

—Volumetric latent heat of fusion, 

L = 144(100)(0.15) = S160 Btu/ft 3 . (eq. 3-2). 

— Average volumetric heat 
capacity, 

C avg - 100[0.17++(0.75 X 0.15)] 

= 28.3 Btu/ft 3 °P. (eq 3-3) 

—Average thermal conductivity, 
K f = 0.80 Btu/ft hr °F (fig. 2-3) 
K u = 0.72 Btu/ft hr °F (fig. 2-4) 
K 
hr 
The % coefficient is as follows: 

—Average surface temperature 
differential, 

v a = nF/t = 2500/160 

= 16.6°F (16.6°F below 32°F). (eq 3-4) 

— Initial temperature differential, 

v_ - MAT -32 » 37.2-32.0 - 



= 1/2 (K u + K f ) = 0.76 Btu/ft 
F. 



5.2°F (5.2° above 32°F). 

— Thermal ratio, 

a - v Q /v s - 5.2/16.6 - 0.33. 

—Fusion parameter, 

fj - v s (C/L) - 16.6(28.2/2160) - 0.20. 

—Lambda coefficient, 



(eq 3-5) 



(eq 3-6) 



(eq3-7) 



A - 0.89 (fig. 3-1). (eq 3-8) 

Estimated depth of frost penetration, 



1 48 K nF 



• 0.89 / 



48(0.76)(2500) 



5.8 ft. 



2160 



(eq 3-9 



3-4. Multilayer soils. 

A multilayer solution to the modified 
Berggren equation is used for non- 
homogeneous soils by determining 
that portion of the surface freezing (or 
thawing) index required to penetrate 
each layer. The sum of the thicknesses 
of all the frozen (or thawed) layers is 
the depth of freeze (or thaw). The 
partial freezing (or thawing) index re- 
quired to penetrate the top layer is 
given by 



L,d 



i"i 



Rt 



F x (or I x ) - 



2 
24 A l 



(eq 3-10) 



3-3 



*TM 5-852-6/AFR 88-19, Volume 6 

where 

d x = thickness of first layer (ft) 
Rj = d/Kj = thermal resistance of 
first layer. 

The partial freezing (or thawing) index 

required to penetrate the second layer 

is 

L 3 d 2 R 3 

P 3 (or I 2 ) - „ (H 1 + — ) . (eq 3-U) 

24 A | 2 

The partial index required to penetrate 
the n m layer is: 



n r 



Fr,(orI n )- 



24A» 



(XR + 



(eq3-12) 



where IR is the total thermal resistance 
above the n th layer and equals 



R, + Ro + Rq ... + R, 



nl 



(eq 3-13) 



The summation of the partial indexes, 



Fj + P 2 + P 3 ... + F n <or I, + I 2 + I 3 . 



(eq 3-14) 



is equal to the surface freezing index 
thawing index). 

a. In this example, determine the 
depth of thaw penetration beneath a 
bituminous concrete pavement for the 
following conditions: 

—Mean annual temperature (MAT) 
- 12°F. 

— Air thawing index (I) 

= 780 degree-days. 
—Average wind speed in summer 

= 7.5 miles per hour (mph). 
—Length of thaw season (t) 

* 105 days. 
—Soil boring log: 



Since a wind speed of 7-1/2 mph results 
in an n-f actor of 2.0 (fig. 2-10), a surface 
thawing index nl of 1560 degree-days 
is used in the computations. The v 8 , v 
and a values are determined in the same 
way as those for the homogeneous case: 

v s » 1560/105 = 14.8°F (eq 3-15) 

v Q - 13.0 - 32.0 » 20.0°F (eq 3-16) 

„ » 20.0/14.8-135. (eq3-17) 

The thermal properties C, K and L of 
the respective layers are obtained from 
figures 2-1 through 2-8. 

b. Table 3-1 facilitates solution of 
the multilayer problem, and in the fol- 
lowing discussion, layer 3 is used to 
illustrate quantitative values. Columns 
9, 10, 12 and 13 are self-explanatory. 
Column 11, L", represents the average 
value of L for a layer and is equal to 
ZLd/Id (2581/5.0 - 517). Column 14, ~C, 
represents the average value of C and 
is obtained from ZCd/Id (145/5.0 = 29). 
Thus TZ and C~ represent weighted 
values to a depth of thaw penetration 
given by Zd, which is the sum of all 
layer thicknesses to that depth. 
The fusion parameter // for each layer 
is determined from 



v_ (C /L) - 14.8 (29/517) - 0.83 . 



(eq 3-18) 



The A coefficient is equal to 0.508 from 
figure 3-1. Column 18, R n , is the ratio 
d/K and for layer 3 equals (3.0/2.0) or 
1.5. Column 19, ZR, represents the sum 
of the R n values above the layer under 
consideration. Column 20, ZR + (R n /2), 
equals the sum of the R n values above 
the layer plus one-half the R n value of 
the layer being considered. Foi ±ay er 3 
this is [1.32 + (1.50/2)] = 2.07. Column 21, 
nl, represents the number of degree- 
days required to thaw the layer being 
considered and is determined from 





Depth 




Dry unit weight 


Water content 


Layer 


(ft) 
0.0-0.4 


Material" 


(lb/ft 3 ) 


(%) 


1 


Asphaltic concrete 


138 





2 


0.4-2.0 


GW-GP 


156 


2.1 


3 


2.0-5.0 


GW-GP 


151 


2.8 


4 


5.0-6.0 


SM 


130 


6.5 


5 


6.0-8.0 


SM-SC 


122 


4.6 


6 


8.0-9.0 


SM 


116 


5.2 



"in accordance with Unified Soil Classification System. 
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Table 3-1 . Multilayer solution of modified Berggren equation. 
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Ld 



R T 



nl 



(IR + ). 

24^2 2 



For layer 3, 



(610K3.0) 



nl 



(2.07) 



(eq3-19) 



(eq 3-20) 



3 r> 

24(0.508) 3 

= 612 degree-days 

The summation of the number of 
degree-days required to thaw layers 1 
through 4 is 1297, leaving (1560 ■ 1297 =) 
263 degree-days to thaw a portion of 
layer 5. A trial and -error method is used 
to determine the thickness of the 
thawed part of layer 5. First, it is as- 
sumed that 1.0 feet of layer 5 is thawed 
(designated as layer 5a). Calculations 
indicate 465 degree-days are needed 
to thaw 1.0 foot of layer 5 or (465 - 263 =) 
202 degree-days more than available. 
A new layer, 5b, is then selected by the 
following proportion 

(263/465)1.0 = 0.57 ft (try 0.6 ft). (eq 3-31) 

This new thickness results in 260 
degree-days required to thaw layer 5b 
or 3 degree-days less than available. 
Further trial-and-error is unwarranted 
and the total estimated thaw penetra- 
tion would be 6.6 feet. A similar tech- 
nique is used to estimate frost penetra- 
tion in a multilayer soil profile. 

3-5. Effect of snow and vegetative cover. 

Thermal properties of snow and vegeta- 
tive covers are extremely variable in 
both time and space. Both materials 
tend to act as insulators and retard 
heat transfer at the air-ground inter- 
face. In f reeze-thaw computations, snow 
and vegetative surface materials are 



equation, with snow cover thickness 
estimated seasonally. The tabulation 
below presents average thermal proper- 
ties of snow applicable for calculation 
in the noted regions if a better data 
base is not available. In the absence of 
site-specific data, figures 2-5 and 2-6 
should be used to estimate the thermal 
conductivities of vegetative surface 
cover. 

3-6. Surface temperature variations. 

The temperatures at the air-ground 
interface are subject to daily and sea- 
sonal fluctuations. Precipitation, insola- 
tion, air temperature and turbulence 
contribute to these variations in sur- 
face temperature. To facilitate mathe- 
matical calculations, two assumptions 
are commonly made regarding the tem- 
peratures at the upper boundary: 1) a 
sudden step change occurs in surface 
temperature or 2) the surface tempera- 
ture change is sinusoidal. The sinu- 
soidal variation of temperature over a 
year closely approximates actual condi- 
tions; however, it is amenable to hand 
calculations only if latent heat effects 
are negligible. Solutions and examples 
for both conditions are given below. 

a. Suddenstep change. This involves 
a sudden change in the surface tem- 
perature of a mass that was initially at 
a constant, uniform temperature. The 
sudden step change was used to estab- 
lish the boundary conditions for heat 
f low in the modified Ber ggren equation 
given in paragraph 3-2. If the influence 
of latent heat is not involved, or is as- 
sumed negligible, the following equa- 
tion may be used: 



eateu as sepfcu'tut; it*y tjx-s in wit? uiuiu- 
y er solution of the mod if ied Ber ggren 


T (x.t) = 


= T s + (T - T s ) erf ( 


2 Vat 


) (eq 3-22) 


Region 


Unit 
weight 
(lb/ft 3 ) 


K 
(Btu/ft) 
hr°F) 

0.11 


C 

(Btu/ft 3 ° F) 




L 
(BtU/ft 3 ) 


Interior Alaska 


16 


8 




2300 


Canadian Archipelago, 
N. Alaskan coast, 
and temperate regions 


20 


0.18 


10 




2880 


Northern Greenland 


22 


0.20 


11 




3170 
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where 
T 



(x,t) 



T c 



T, 



erf 



= temperature at depth x, at 
time (°P) 

= suddenly applied constant 
surface temperature (°F) 

= initial uniform tempera- 
ture of the mass (°P) 

= mathematical expression, 
termed the error function, 
which is frequently used 
in heat flow computations 
(dimensionless) 

= depth below surface (ft) 



a - thermal diffusivity of the 

mass (ft 2 /day = K/C) 
t = time after application of 

sudden change in surface 
temperature (days). 
Figure 3-2 give s the relatio nshi p be- 
tween (x/2 M at)anderf(x/2 V~~at).The 
expression for the error function is 
shown in appendix A and figure 3-2. 
In this example of a sudden step change, 
a highly frost-susceptible subgrade is 
covered with a 2-foot thick, non-frost- 
susceptible gravel pad. Both soils are 
at an initial temperature of 20°F. If the 
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Figure 3-2. Relationship between (x/2 Y"at) and erf (x/2 Vat). 
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surface of the gravel is suddenly heated 
to and maintained at 70°P for a number 
of days, estimate the temperature at 
the gravel-subgrade interface after one 
day. The gravel material is very dry 
and latent heat may be ignored. The 
thermal conductivity of the gravel is 
1.0 Btu/f t hr °F and the volumetric heat 
capacity is 25 Btu/f t 3 . The thermal dif- 
fusivity of the gravel is (K/C - 1.0/25) 
== 0.04 ft 2 /hr., o r 0.96 ft 2 /day, and 
x/2 Vat = (2.0/2 V096X 1) = 1.02. From 
figure 3-2, erf (x/2 Vat) is equal to 0.85, 
and the interface temperature T is 
[70 + (20 - 70)0.85] = 27.5°F. 

b. Sinusoidal change. A surface tem- 
perature variation that is nearly sinu- 
soidal repeats itself periodically for a 
surface exposed to the atmosphere. For 
most problems in this manual, the sinu- 
soidal variation of concern occurs over 
an annual cycle. If latent heat is not 
involved or is assumed negligible, the 
following equation may be used: 



A x = A Q exp (-x 



" aP 



(eq 3-23) 



where 

A= amplitude of temperature wave 
at depth x (°F). 

A= amplitude of the surface tem- 
perature wave above or below 
the average annual tempera- 
ture (°F) 

x = depth below surface (ft) 

a as thermal dif fusivity of the mass 
(ft 3 /day) 

P = period of sine wave (365 days). 
The sinusoidal temperature pattern is 
assumed to exist at all levels to a depth 
where there is no temperature change. 
The temperature waves lag behind the 
surface wave, and the amplitude of the 
sinusoidal waves decreases with depth 
below the surface. The pha se lag i s 
determined by t x = (x/2) ( "V 365/ „-a). 
Typical temperature-time curves for a 
surface and at a depth x are shown in 
figure 3-3. 

In the following example of a sinu- 
soidal temperature change, the surface 
temperature of an 8-foot- thick concrete 
slab varies from 60° to -40°F during the 
year. Determine the maximum tem- 
perature at the base of the slab as- 
suming a diffusivity of 1.0 ft 2 /day for 
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the concrete. The average annual tem- 
perature is [60 + (-40)] /2 = 10°F and the 
surface amplitude is (60 - 10) - 50°F. 
The amplitude at an 8-foot depth equals 



1 50 exp [-8 



4 



, 5Oe .0.742 



(1.0)(365) 

- 24°F. (eq 3-24) 

The maximum temperature at 8 feet is 
(10 + 24) = 34°F. The time lag t x between 
the maximum temperature at the sur- 
face and 8 feet is 



^ \ 



365 



=43 days 



ff (1.0) (about 6 weeks). (eq3-25) 

(Note: latent heat would increase the 
time lag and decrease the amplitude.) 

c. Nonuniform layers. The method 
of equivalent thickness is used to find 
the temperature at a point below a 
number of layers of different thermal 
properties. This technique assumes a 
negligible effect of latent heat, and in- 
volves reduction of each layer to an 
equivalent material thickness by set- 
ting the ratio of the thicknesses equal 
to the ratio of square roots of the ther- 
mal diffusivities. For example, deter- 
mine the equivalent gravel thickness 
for the three layers shown, assuming 
all materials are unfrozen. 
The following table shows that 4.75 
feet of the nonuniform materials can 
be considered equivalent to 5.4 feet of 
gravel for heat- flow purposes. This 
equivalent thickness and the thermal 
diffusivity of the gravel are used to 
calculate temperatures at the base of 
the gravel layer by either the step- 
change or sinusoidal method. 

3-7. Converting indexes into equivalent sine 
wave of temperature. 

Some problems may require the use of 
the sinusoidal temperature variation 
technique, given only the freezing or 
thawing indexes and the average 
annual temperature. These indexes 
may be converted into a sine curve of 
temperature to give the same index 
values and the sam emean temperature. 
For example, convert the monthly aver- 
age temperature data for Fairbanks, 
Alaska, shown in figure 2-9 into an 
equivalent sine wave. The relationship 
between the sinusoidal amplitude, 
freezing index, thawing index and aver- 
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SURFACE TEMPERATURE 




FREEZING INDEX AT DEPTH (X) 



TIME 

(U.S. Army Corps of Engineers) 

Figure 3-3. Sinusoidal temperature pattern. 



Dry Unit 
weight, 

Yd 

Material (lb/ft 3 ) 


Water 

content, 

w (%) 

2 

4 


Thermi 
conducts 

K 
(Btu/ft hr 

1.0 

0.8' 

1.5" 


»/ 
<ity, 

°F) 


Volumetric 

heat 
capacity, 

C 
( Btu/ft 3 °F) 


Thermal 

diffusivity, 

a = K/C 

(ft 2 /hr) 


Concrete 

Sand 120 

Gravel 135 




33.0 

23 1 
28 fx 


0.033 
0.035 
0.054 


"From fiqure 2-2. 

tC = y d (0.17 + w/100). 

Material 


Thickness 
(ft) 




Thermal 

diffusivity 

(fi 2 /hr) 


V^rr? 


Equivalent gravel 
thickness (ft) 


Concrete 

Sand 

Gravel 


1.75 
0.50 
2.50 




0.033 
0.035 
0.054 


1.3 
1.2 
1.00 


2.3(1.3X1.75) 
0.6 
2.50 


Total thickness 


4.75 




5.4 



The subscript g refers to the gravel layer and the subscript m refers to the other material layer. 
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age annual temperature Is shown in 
figure 3-4. 

The average temperature at Fairbanks 
for October 1949 to September 1950 
was 27°F. By use of the freezing index 



of 5240 degree-days, the sinusoidal 
amplitude is found to be 37.0°F. The 
equation of the sine wave is 

T = 27.0 - 37.0 sin Zwtt (eq 3-26) 

= 27.0 - 37.0 sin 0.0172 t (radians) 
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Figure 3-4. Indexes and equivalent sinusoidal temperature. 
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where 

f = frequency, 1/365 cycles per day 
t = time from origin in days. (Origin 
of curve is located at a point 
where T intersects the average 
annual temperature on its way 
downward toward the yearly 
minimum.) 

If the -thawing index of 3240 degree- 
days had been used, the sinusoidal tem- 
perature amplitude would be 35.5°P. 
The actual temperature curve for Pair- 
banks, Alaska, and the equivalent sine 
wave computed from the freezing index 
are plotted in figure 3-5. This illustra- 
tion makes use of air indexes butasine 
wave could be determined for surface 
indexes by multiplying the air indexes 
with appropriate n-factors. Note that 
the mean annual ground surface tem- 
perature may be substantially different 
(frequently higher) from the mean an- 
nual air temperature because the 
freezing n-f actor is generally not equal 
to the thawing n-factor. If the long- 
term mean monthly temperature had 
been used instead of the average 
monthly temperatures for the 1949-1950 
period, the correlation between the 
actual temperature curve and the equi- 
valent sine curve would practically coin- 
cide, as shown in figure 3-6. 




(U.S. Army Corps of Engineers) 

Figure 3-5. Average monthly temperatures for 

1949-1950 and equivalent sine wave, Fairbanks, 

Alaska. 



3-8. Penetration of freeze or thaw beneath 
buildings. 

The penetration of freeze or thaw 
beneath buildings depends largely on 
the presence or absence of an airspace 
between the building floor and the 
ground as discussed below and in TM 
5-852-4/ APM 88-19, Chapter 4. 

a. Building floor placed on ground. 
When the floor of a heated building is 
placed directly on frozen ground, the 
depth of thaw is determined by the 
same method as that used to solve a 
multilayer problem when the surface 
is exposed to the atmosphere, except 
that the thawing index is replaced by 
the product of the time and the dif- 
ferential between the building floor 
temperature and 32°F. For example, 
estimate the depth of thaw after 1 year 
for a building floor consisting of 8 
inches of concrete, 4 inches of insula- 
tion and 6 inches of concrete, placed 
directly on a 5-foot- thick sand pad over- 
lying permanently frozen silt for the 
following conditions: 

—Mean annual temperature = 20°P. 

— Building floor temperature = 65°F. 

—Sand pad: y d - 133 lb/ft 3 , w ■ 5%. 

—Frozen silt: y d ■ 75 lb/ft 3 , w = 45%. 

—Concrete: K = 1.0 Btu/ft hr 
°P, C = 30 Btu/ft 3 °F. 

—Insulation: K = 0.033 Btu/ft hr °F, 
C = 1.5 Btu/ft 3 °F. 
The resistances of the three floor 
layers are in series, and the floor re- 
sistance R f is the sum of the three layer 
resistances: 



Ft* 



8 



d 
k 
11.2 ft* ft ° 



, (12)(1.0) 



(12)(0.033) 



(12)(1.0) 
(eq 3-27) 



The average volumetric heat capacity 
of the floor system is 

„ (30)(8) + (jLS)H) + (30X6) 

1 8+4+6 

= 23.7 Btu/ft 3 °F (eq 3-28) 

The solution to this problem, shown in 
table 3-2, predicts a total thaw depth of 
7.8 feet. This solution did not consider 
edge effects, i.e.,a long narrow building 
will have lesser depth of thaw than a 
square building with the same floor 
because of the difference in lateral heat 
flow. 
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NOTE: Areas between curves and 32'F line 
represent the indexes. 



SEPT I OCT | NOV 1 DEC 1 JAN | FEB | MAR | APR | MAY j JUNE | JULY [ AUS | SEPT | OCT 



(U.S. Army Corps of Engineers) 

Figure 3-6. Long-term mean monthly temperatures and equivalent sine wave, Fairbanks, Alaska. 



b. Airspace below building 
(1) An unskirted airspace between 
the heated floor of building and the 
ground will help prevent degradation 
ofunderlyingpermafrost-Theairspace 
insulates the building floor from the 
ground and acts as a convective pas- 
sage for flow of cold air that dissipates 
heat from the floor system and the 
ground. The depth of thaw is calculated 
by means of the modified Berggren 
equation for either a homogeneous or 
multilayer ed soil system, as applicable. 
An n-factor of 1.0 is recommended to 
determine the surface thawing index 
beneath the shaded area of an elevated 
building. 
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(2) There is no simple mathemati- 
cal expression for analyzing the heat 
flow in a ventilated floor system that 
has ducts or pipes installed within the 
floor or at some depth beneath the floor, 
with air circulation induced by stack 
effect. The depth to which freezing 
temperatures will penetrate is com- 
puted with the modified Berggren equa- 
tion, except that the air freezing index 
at the outlet governs. This index is 
influenced by a number of design vari- 
ables, i.e., average dally air tempera- 
tures, inside building temperatures, 
floor and duct or pipe system design, 
temperature and velocity of air in the 
system, and stack height. Cold air 
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Table 3-2. Thaw penetration beneath a slab-on-grade building constructed on permafrost. 
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passing through the ducts acquires 
heat from the duct walls and experi- 
ences a temperature rise as it moves 
through the duct, and the air freezing 
index is reduced at the outlet. Field 
observations indicate that the inlet air 
freezing index closely approximates 
the site air freezing index. The freezing 
index at the outlet must be sufficient 
to counteract the thawing index and 
ensure freeze-back of foundation soils. 

(3) As an example, determine the 
required thickness of a gravel pad be- 
neath the floor section shown in figure 
3-7 to contain all thaw penetration. Also 
determine the required stack height to 
ensure freeze-back of the pad on the 
outlet side of the ducts. The conditions 
for this example follow: 

— Duct length, 1 - 220 ft. 

— Gravel pad: y d - 125 lb/ft 3 , 
w = 2.5%. 



— Outlet mean annual tempera- 
ture - 32°F 0, w - 2.5%. 
(conservative assumption). 

— Minimum site freezing index 

- 4000 degree-days. 

— Freezing season = 215 days. 

— Thawing season = 150 days 
(period during which ducts are 
closed). 

— Building floor temperature 

- 60°F. 

— Thermal conductivity of con- 
crete, K c ■ 1.0 Btu/f t hr °P. 

— Thermal conductivity of insula- 
tion, Kj ■ 0.033 Btu/f t hr °° 00 F. 

(a) The required thickness is 
determined by the following equation, 
derived from the modified Berggren 
equation: 



48t*I|, 



KR, 



1 + 



KL(Rf)'' 



1] 



(eq 3-29) 



Figure 3-7. Schematic of ducted foundation. 
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where 

K = average thermal conductivity of 
gravel 
- 1/2(0.7 + 1.0) = 0.85 Btu/f t hr °F 
R T = thermal resistance of floor 



system 

18 



12 



12(1.0) 12(0.033) 12(1.0) 

= 12.5 ft 2 hr °F/Btu (eq 3-30 

(In the computations the dead airspace 
is assumed equivalent to the thermal 
resistance of concrete of the same 
thickness.) 

X = factor in modified Berggren 
equation = 0.97 (conservative 
assumption) 
I r = thawing index at floor surface 
= (60 - 32)(150) = 4200 degree- 
days 
L=latent heat of gravel = 
144(125X0. 025) » 450 Btu/ ft3 

then / (48)(0.97) 2 (4200) 

X - (0.85)(12.5) [ /l + 1] 

\ (0.85)(460)(12.5) 2 

= H.O ft. (eq 3-31) 

(b) Thus the total amount of heat 
to be removed from the gravel pad by 
cold-air ventilation during the freezing 
season with ducts open is equal to the 
latent and sensible heat contained in 
the thawed pad. The heat content per 
square foot of pad is determined as 
follows: 

— Latent heat, (X)(L) - (H.O)(450) 
- 4950 Btu/ft 2 

— Sensible heat (10 percent of 
latent heat, based upon ex- 
perience) = 495 

— Total heat content: 

5445 Btu/ft 2 
The ducts will be open during the 
freezing season (215 days), and the aver- 
age rate of heat flow from the gravel 
during this season is equal to 5445/215 
X 24 = 1.0 Btu/ft 2 hr. The average thaw- 
ing index at the surface of the pad is 



LX* 



<450)<U.Or 



48A 2 K 48(0.97) 2 (0.85) 



= 1420 degree-days. 



(eq 3-32) 



This thawing index must be com- 
pensated for by an equal freezing index 
at the duct outlet on the surface of the 
pad to assure freeze-back. The average 



pad surface temperature at the outlet 
end equals the ratio 

Required freezing index 

Length of freezing season 
1420 



215 
- 6.6°F below 32°F or 25.4°F. 

The inlet air during the freezing season 
has an average temperature of 

Air freezing index 4000 



Length of freezing season 
■ 18.6°F below 32°F or 13.4°F. 



215 



Therefore, the average permissible tem- 
perature rise T R along the duct is 
(25.4 - 13.4) = 12.0°F. 

(c) The heat flowing from the 
floor surface to the duct air during the 
winter is equal to the temperature dif- 
ference between the floor and duct air 
divided by the thermal resistance be- 
tween them. The thermal resistance R 
is calculated as follows: 



x^ 



x, 



R - + +■ 

K„ K, 



14 



(12X1.0) 



= 12.3 hr ft* 5 "F/Btuteq 3-33) 



(12X0.033) 1.0 

where 

X c = thickness of concrete (ft) 
X t = thickness of insulation (ft) 
h rc = surface transfer coefficient 
between duct wall and duct 
air 
(For practical design, h rc = 1.0 Btu/ft 2 
hr °F and represents the combined 
effect of convection and radiation. At 
much higher air velocities, this value 
will be slightly larger; however, using 
a value of 1.0 will lead to conservative 
designs). The average heat flow be- 
tween the floor and inlet duct air is 
[(60 - 13.4)/12.3] - 3.8 Btu/ft 2 hr, and be- 
tween the floor and outlet duct air is 
[(60 - 25.4)/12.3] - 2.8 Btu/ft 2 hr. Thus 
. the average rate of heat f low from the 
gravel pad to the duct air is l.O Btu/ft 2 
hr. The total heat flow to the duct air 
from the floor and gravel pad is 
(3.3 + 1.0) = 4.3 Btu/ft 2 hr. The heat flow 
to the duct air must equal the heat 
removed by the duct air* 

Heat added ■ heat removed 



#|m-60VA d /,c p T R . 



(eq 3-34) 
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Thus the average duct air velocity 
required to extract this quantity of heat 
(4.3 Btu/ft 2 hr) is determined by the 
equation: 

v = ft/minute (eq3-35) 

60A d pC p T R 

wher« 

= total heat flow to duct air (4.3 

Btu/ft 2 hr) 
i - length of duct (220 ft) 
m = duct spacing (2.66 ft) 
A d = cross-sectional area of duct 

(1.58 ft 2 ) 
<p -density of air (0.083 lb/ft 3 

[figure 3-10]) 
c p = specific heat of air at constant 

pressure (0.24 Btu/lb °F) 
T R - temperature rise in duct air 

(12°F). 



Substitution of appropriate values 
gives a required air velocity 

<4.3)<220)<2.66) 

V= 

(G0)(1.58)(0.083)<0.24)(12.0) 
■ 111 ft/minute. <eq 3-36) 

(d) The required air flow is ob- 
tained by a stack or chimney effect, 
which is related to the stack height. 
The stack height is determined by the 
equation 

h d -h v + h f (eq3-37) 

where 

p e H(T C - T ) 

h rt ■ inches of water 

5.2(T c + 460) 
(natural draft head) 

p ■ density of air at average duct 

temperature (lb/ft 3 ) 
e * efficiency of stack system (%). 

This factor provides for 



Figure 3-8. Properties of dry air at atmospheric pressure. 
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friction losses within the 

chimney 
H = stack height (ft) 
T c = temperature of air in stack 

(°F) 
T = temperature of air 

surrounding stack (°F) 

V 2 ^ 



h v =( 



) 2 inches of water 

(velocity head) 



4000 
V = velocity of duct air (ft/minute) 

h f = f — h v inches of water 
D e (friction head) 
f = friction factor 

(dimensionless) 
£ e = equivalent duct length (ft) 
D e = equivalent duct diameter (ft). 
The technique used to calculate the 
friction head is 



4(cross-sectional area 
of duct in ft 3 ) 



4(1.58) 



*V 



2 t 18 + 30 ._, 
— ;< 5~ +12) 



perimeter of duct in ft 

12 * 2 

= 1.22 ft. (eq 3-38) 

The equivalent length of the duct is 
equal to the actual length i s plus an 
allowance I b for bends and entry and 
exit. Each right-angle bend has the 
effect of adding approximately 65 
diameters to the length of the duct, and 
entry and exit effects add about 10 
diameters for each entry or exit. In this 
example the total allowance i b for these 
effects is [2(65 + 10) =] 150 diameters, 
which is added to the length of the 
straight duct. The estimated length of 
straight duct^ s is 

5 ft (assumed inlet open length) 
220 ft (length of duct beneath 

floor) 
15 ft (assumed stack height) 



240 ft 

j _ » 240 + (150 X 1.22) - 423 ft. 



(eq 3-39) 



The friction factor f is a function of 
Reynolds number N R and the ratio 
e/D e . A reasonable absolute roughness 
factor e of the concrete duct surface is 
0.001 feet, based on field observations. 
Suggested values of e for other types 
of surfaces are given in the ASHRAE 
Data and Guide Book. The effect of 
minor variations in e on the friction 



factor is small, as noted by examining 
the -equation below used to calculate 
the friction factor f . Reynolds number 
is obtained from the equation 

V(a' + 0.25 D e ) 



N 



R 



(eq 3-40) 



where 



(111 X 60)(1.0 + 0.25 X 1.22) 



N R - 



■ 17,700 



0.49 

V » average duct velocity (ft/hr) 
a ■ shortest dimension (ft) 
V » kinematic viscosity 

(ft 2 /hr at 19.4°F [fig. 3-8]). 

The friction factor f* is obtained by 
solving the equation 



F - 0.0055(1 + (20.000 X + 

D e 
0.001 

- 0.0055[1 + (20.000 X 

1.22 



io c 



H 



R 



) 1/3 1 



10' 



,8 



a/3i 



17.700 
- 0.0285. (eq 3-41) 

Therefore, the friction head is 

'a 

(eq3-42) 



\±/Ol 



h.» f x -r— 1 x h, 
' D e " 



- 0.0285 X 



423 
1.22 



Xh„»Mh, 



The draft head required to provide the 
desired velocity head and to overcome 
the friction head is furnished by the 
chimney or stack effect- The draft head 
h d is obtained as follows: 



h„ + 9.8h v 



(eq 3-43) 



10.8 



10.8 (- 



* 10.8 (- 



"v 

V 

4000 
HI 

4000 



*r 



8.31X10" 3 inches of water. 



The stack height required to produce 
this draft head is 



H 



5.2 h d (T c + 460) 



p £ (TVT„) 



(eq 3-44) 



(5.2)(8.31X10 3 )(25.4 + 460) 

(0.083)(0.80)(25.4 - 13.4) 
- 26 ft 

where 

p m 0.083 lb/ft 3 
T c »25.4°F 
T - 13.4°F 

e « 80% (found to be a reasonable 
design value based on observa- 
tions over an entire season) 
h d - 8.31X10- 3 inches of water. 
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(e) If the stack is too high for the 
struc ture, a greater thickness of insula- 
tion could be used. In this example, the 
effect of increasing the insulation thick ■ 
ness by one-half would result in lower- 
ing the stack height by five-eighths. 

(f) This first approximated stack 
height is next incorporated in the cal- 
culation of the length of straight duct 
i s , and the newly obtained t e is used to 
recalculate the friction head h P By trial - 
and -error, the final calculated stack 
height is found to be 26.5 ft. 

(g) The stack height is an im- 
portant variable because an increase 
in stack height will increase the duct 
airflow. Circulation of air through the 
ducts results from 1) a density dif- 
ference between the air inside the duct 
and that outside the building, 2) apres- 
sure reduction at the outletend attrib- 
utable to the stack effect, 3) a positive 
pressure head at the inlet end when 
wind blows directly into the intake 
stack opening, and 4) a negative pres- 
sure head at the outlet when wind 
passes over the exhaust stack opening. 
Draft caused by wind is highly erratic 
and unpredictable and should not be 
considered in design; however, the 
vents should be cowled to take ad- 
vantage of any available velocity head 
provided by the wind. If sufficient air 
cannot be drawn through ducts by 
natural draft, mechanical blowers 
could be specified or consideration 
given to alternating air flow in the ducts. 

3-9. Use of thermal insulating materials. 
An insulating layer may be used in 
conjunction with a non-frost-suscep- 
tible material to reduce the thickness 
of fill required to keep freezing or 
thawing temperatures from pene- 
trating into an underlying frost-suscep- 
tible soil. As in the example of para- 
graph 3-8a, the thermal resistance of 
the pavementand insulation layers are 
added to obtain total resistance, and 
the latent heat effect of a combined 
pavement and insulation layer is as- 
sumed negligible. (TM 5-818-2/ AFM 88- 
6, Chap. 4 discusses in detail the design 
of insulated pavements.) If the insu- 
lating material will absorb water, its 
insulating effectiveness will be re- 
duced considerably (as discussed in 



TM 5-852-4/ AFM 88-19, Chap. 4). Limited 
field tests indicate that the heat-flow 
resistance of a portland-cement- 
concrete pavement overlying a high- 
quality insulating layer is more compli- 
cated than simple addition of resist- 
ances, but until sufficient data are 
obtained for validation, treatment of 
resistances in series is recommended. 
a. Example. A pavement consists of 
14 inches of portland-cement concrete 
placed on a 6-foot gravel base course. 
Frostpenetrated3.2 feet into the under- 
lying silt subgrade. Determine the thick- 
ness of insulation required to prevent 
frost penetration into the subgrade for 
the following conditions. 

— Mean annual temperature = 
35.3°F. 

— Air freezing index = 
3670 degree-days 

— Freezing season = 170 days 

— Concrete: K ■ 1.0 Btu/ft hr °F, 
C - 0.30 Btu/ft 3 °F 

— Insulation: K = 0.024 Btu/ft hr 
°F, C = 0.28 Btu/ft 3 °F 

— Gravel base: y d = 130 lb/ft 3 , 
w = 4% 

— Silt subgrade: y d = 100 lb/ft 3 , 
w = 10%. 

Surface freezing index = 0.75 X 3670 = 
2752 degree-days. From the known 
data 



2752 



170 



16.2°F 



35.3 - 32.0 = 3.3F 
3.3 



16.2 



= 0.20 



(eq 3-45) 



(eq 3-46) 



(eq 3-47) 



b. Trial i.Usea2-inchlayerofinsula- 
tion and a 6-inch concrete leveling 
course. 

— Pavement section: 14 Inches of concrete 
2 inches of insulation 
6 inches of concrete 
leveling course 





22 inches total 


R p 


d 14 2 6 


K 12 X 1.0 12 x 0.024 12 X 1.0 


= 


8.60 hr ft 2 °F/Btu (eq 3-48) 




(14X30) + (2X0.28) + (6X30) 


C P 


22 




= 27.3 Btu/ft 3 °F (eq 3-49) 
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The calculation appears in table 3-3 
and indicates that this pavement sec- 
tion has an excess of (2752 - 3033 =)481 
degree-days to prevent frost pene- 
tration into the silt subgrade. 

c. Trial 2. Use a 1.5- inch layer of 
insulation and a 6-inch concrete level- 
ing course. 

—Pavement section: 14 inches of concrete 
1.5 inches of insulation 
6 Inches of concrete 
leveling course 



14 



1.5 



R P = 



12 X 1.0 
(eq 3-50) 



21.5 inches total 



K 12X1.0 12X0.024 

6.86 hr ft 2 °F/Btu 

(14X30) + (1.5X0.28) + (6X30) 

C p = 

21.5 
■ 27.9 Btu/ft 3 °F (eq3-51) 

The calculation (see table 3-4) indicates 
that this pavement section will not 
prevent frost penetration into the silt 
subgradeas(2752 - 2553 =) 199 degree- 
day s remain for subgrade penetration. 
The 2-inch thickness of insulation is 
therefore required. 
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Table 3-3. Insulated pavement design, no frost penetration. Table 3-4. Insulated pavement design, frost penetration. 
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CHAPTER 4 
TWO-DIMENSIONAL RADIAL HEAT FLOW 



4-1. General. 

Radial flow of heat Is considered in 
thermal problems involving the design 
of pile foundations in permafrost (TM 
5-852-4/ AFM 88-19, Chapter 4), the con- 
struction of utility supply lines for the 
transport of water and sewage in perma- 
frost areas and seasonal frost areas 
(TM 5-852-5/ AFR 88-19, Volume 5), and 
the design of artificially frozen ground 
for retaining structures during con- 
struction. A number of the basic con- 
cepts and techniques used to calculate 
radial heat flow from cylindrical sur- 
faces are discussed below. 

a. Thermal Resistance. In analyzing 
heat flow for areas with cylindrical 
cross sections, the effective thickness 
for radial flow f rom aunit length of the 
cylinder is 



Zn 



In 



or 0.367 log 



where 



rj = inside wall radius (ft) 
r 3 = outside wall radius (ft). 
The thermal resistance R is equal to 
the effective thickness divided by the 
conductivity of the material between 
the two radii. As an example, a concrete 
conduit with a wall thickness of 6 
inches and an inside diameter of 10 feet 
is surrounded by 4 inches of cellular 
glass insulation and 4 feet of dry gravel. 
Calculate the thermal resistance be- 
tween the inside concrete wall and the 
outer edge of the gravel material. The 
following thermal conductivities are 
given: 

— Concrete, K = 1.00 Btu/f t hr °F. 

— Insulation, K = 0.033 Btu/ft hr 

°F. 

— Gravel, K = 1.5 Btu/ft hr °F. 

Let (see fig. 4-1 for values of r x - r 4 ) 
r t = radius to inner wall of conduit 
r = radius to outer wall of concrete 



r 3 = radius to outer edge of insula- 
tion 

r 4 = radius to outer edge of gravel 
K x 3 = thermal conductivity of concrete 
K 2-3 = thermal conductivity of insula- 
tion 

K 3 . 4 - thermal conductivity of gravel. 



R = 0.367 ( 



*l-2 



log 



K, 



log 



(eq 4-1) 



2-3 



K, 



3-4 



log 



5.5 



1 5.83 1 9.83 

= 0.367 ( — log — ■+ log — + — io« — ) 

1.0 5.0 0.033 U 5.5 1.5 S 5.83 ' 

-0.367 (0.041 + 0.767 + 0.151) =0.352 ft. 2 hr °F/Btu. 



If the temperature in the conduit were 
45°F and the temperature at the outer 
face of the gravel were 35°F, the heat 
flow per linear foot of conduit would 
equal 



0.352 



(45 - 35) = 28.4 Btu/hr. 



(eq 4-2) 



b. Temperature field surrounding 
a cylinder. The sudden or step change 
in surface temperature discussed for 
semi-infinite slabs in paragraph 3-6a 
has application to heat- flow problems 
associated with pile foudations in 
permafrost. A mathematical solution 
is available for the problem where the 
surface temperature of a cylinder is 
suddenly changed from the uniform 
temperature of the surrounding 
medium, as long as there is no phase 
change. Figure 4-2 is used to determine 
the temperature T at a distance r from 
the center of a cylinder of radius rj, at a 
time t after the surface temperature of 
the cylinder is changed from T Q to T s . 
The temperature T Q represents the uni- 
form temperature of the medium prior 
to the sudden change in surface 
temperature. . 
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Figure 4-1. Illustration for example in paragraph 4- 1 a. 



4-2. Pile installation in permafrost. 

At many arctic and subarctic sites, pile 
foundations are commonly placed in 
preaugered holes, and the annular 
space between the oversized hole and 
pile is backfilled with a slurry of soil 
and water. The tangential adfreeze 
strength of the pile is principally a 
function of the bond developed be- 
tween the pile and the frozen slurry. 
Dissipation of the sensible and latent 
heat of the slurry into permafrost is a 
major design factor because construc- 
tion scheduling depends upon the time 
required for slurry freeze-back. Pile 
spacing is important as each pile adds 
heat, i.e., piles spaced too closely may 
increase permafrost temperatures, 
with a reduction of pile adfreeze 
strength and an increase in freeze- 
back time. This is particularly true in 
relatively warm permafrost (above 
30°F). The following discussion as- 
sumes first that the slurry will freeze 
back naturally because of heat transfer 
between the surrounding permafrost 
and the slurry, second that the time 
required for freeze-back ata particular 
depth is predominantly influenced by 
the permafrost temperature at that 
depth, and third that the permafrost 
does not thaw. Under certain condi- 



tions, artificial refrigeration may be 
required to ensure freeze-back within 
a reasonable time. The volumetric 
specific heat of the slurry and the effect 
of vertical heat flow are assumed to 
have a negligible effect in computing 
required freeze-back time. With proper 
pile spacing, the slurry temperature 
reaches that of the surrounding perma- 
frost in time. Surrounding tempera- 
tures, natural freeze-back time, proper 
pile spacing, and heat removal by refrig- 
eration, as well as heat transfer by ther- 
mal piles, are discussed below. 

a. Surrounding temperatures. The 
increase in permafrost temperatures 
during slurry freeze-back is deter- 
mined using the technique described 
in paragraph 4-lb.. For example, a 
preaugered hole for a pile installation 
is 16 inches in diameter and the slurry 
is placed at a temperature of 32°F. The 
surrounding permafrost has a thermal 
diffusivity of 0.06 ft 2 /hr and an initial 
temperature of 28°F. Calculate the 
ground temperature at a distance of 3 
feet from the center of the pile after an 
elapsed time of 48 hours. Given 



r t = 0.667 ft 

r - 3.00 ft 

a - 0.06 ft 3 /hr 



t - 48 hr 
T = 28°F 
T s = 32°F 
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Figure 4-2, Temperature around a cylinder having received a step change in temperature. 



therefore, 



r l 


0.667 


vir 


"V (0.06)(48) 


r 


3.00 
= = 4.50 



From figure 4-2, 

T-T„ 



0.39 



0.667 



(eq 4.3) 



(eq 4-4) 



: =0.18 

T ~ 0.18 (T s - T )+T Q = 0.18 <32-28)+28 
- a8.7°F. 



(eq4-5) 
(eq 4-6) 
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This technique is also used to predict 
the increase in permafrost temperature 
during slurry freeze-back. Since it as- 
sumes a constant surface temperature 
for a cylinder, it is applicable only to 
the time of freeze-back. After the slurry 
has frozen, the permafrost temperature 
decreases and the model is not valid. 

b. Natural freeze-back time. This 
heat transfer problem assumes a slur- 
ried pile to be a finite heat source inside 



a semi-infinite medium, with asuddenly 
applied constant temperature source 
(32°P) that dissipates radially into 
frozen ground of known initial tempera- 
ture. The general solution for the 
natural freeze-back problem, based 
upon the latent heat content of the 
slurry, is shown in figure 4-3. 

(1) For example, calculate the time 
required to freeze back a 12.5-inch 
diameter timber pile placed in an 18- 



Figure 4-3. General solution of slurry freeze-back. 
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inch hole preaugered in permafrost 
and backfilled with a slurry for the 
following conditions: 

— Permafrost: Sllty sand 

Initial temperature - 27°F 
Dry unit weight - 94 lb/ft 3 
Water content » 25%. 

— Slurry backfill: Silt, water 

Placement tempera- 
ture - 33.5°F 

Dry unit weight - 72 lb/rt 3 
Water content ■ 45%. 

In figure 2-3, the thermal conductivity 
of the permafrost is determined to be 
1.1 Btu/f t hr°F. The volumetric heat eapa- 
city is calculated to be 



25 o 

[94 (0.17 + 0.5 — )] - 27.7 Btu/ft ,J °F 

and the thermal diffusivity to be 

1.1 



(eq 47) 



27.7 



■ 0.0397 ft a /hr. 



(eq 4-8) 



The volumetric latent heat of the slurry 
is 

(144 X 72 X 0.45) - 4670 Btu/ft 3 °F (eq 4-9) 

and the latent heat per linear foot of 
backfill is 



[ — (1.5 3 - 1.04 3 )] 4670 - 4280 Btu. 



(eq 440) 



When figure 4-3 is entered with a value 

of 
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be increased by 50 percent to permit 
an element of safety in the design. 
(Note: the sensible heat introduced by 
the pile and slurry was negligible in 
comparison to the latent heat intro- 
duced by the slurry and was not con- 
sidered in calculations.) 

(2) Permafrost temperature varia- 
tions with depth, as discussed in TM 
5-852-4/ AFM 88-19, Chap. 4, should be 
considered in calculating freeze-back 
time. Figure 4-4 illustrates the effect 
of permafrost temperature on freeze- 
back for the above example. Since heat 
input is governed principally by the 
latent heat of slurry backfill, which is a 
function of slurry volume, moisture 
content and dry unit weight, a family 
of curves relating volumetric latent 
heat of slurry, permafrost tempera- 
tures and freeze-back time may be de- 
veloped for a specific site to account 



Q 



4280 



r 2At 



<27.7)(0.75 3 )<5) 



] -55 



then, 



at 



12.4. 



(eq 441) 



(eq 442) 



The time required to freeze back the 
slurry backfill is 

12.4 X 0.75 2 

■ 176 hours or about 

0.0397 7.3 days. (eq4 13) 

At this time the slurry temperature is 
32°F. Subsequent to freeze-back, the 
temperature of the slurry will continue 
to decrease and will approach the 
permafrost temperature. Ninety per- 
cent of the temperature difference will 
disappear in about twice the time re- 
quired for freeze-back to 32°F. In this 
example, after a period of [7.3 + (2 
X 7.3) =] 22 days, the slurry temperature 
would be approximately [32 - 0.90 (32 - 
27) =] 27.5°F. The time (22 days) should 



Figure 4-4. Specific solution of slurry freeze-back 
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for varying pile shapes and preauger ed 
hole diameters. To minimize the heat 
introduced by the slurry, the water 
content should be the minimum re- 
quired for complete saturation. This 
can be bestaccomplished by backfilling 
with the highest dry unit weightmate- 
rial that can be processed and placed, 
i.e., a well-graded concrete sand with a 
6-inch slump. 

c. Pile Spacing. The effect of pile 
spacing on the overall rise of perma- 
frost temperature resulting from instal- 
lation of piles in preaugered holes is 
found by equating the latent heat of 
slurry backfill with the allowable sen- 
sible heat (temperature) rise of the sur- 
rounding permafrost. For example, 
calculate the minimum allowable pile 
spacing in the preceding example so 
that the permafrost temperature will 
notrise above31°F.The following equa- 
tion, equating the latent heat of the 
slurry to the change of sensible heat in 
a permafrost prism of side S, is used to 
determine the pile spacing: 



V 



(7rr§) + 



(eq 4-14) 



C A T 

where 

S - grid pile spacing (ft) 
r 2 = radius of augered hole (ft) 
Q = latent heat of slurry per lineal foot 
(Btu/ft) 

C = volu me trie heat capac ity of perma- 
frost (Btu/ft 3 °F) 

AT = temperature rise of permafrost 

(°F). 

Substitution of appropriate values from 
the above example and a maximum allow- 
able permafrost temperature rise A T 
of 4°F give a minimum pile spacing S of 



V 



(3.14)(0.75) 2 - 



4280 



5.4 ft. 



(eq 4-15) 



(27.7)(4) 

This spacing may not keep local tem- 
perature from rising to more than 31°F; 
however, it will keep the entire mass of 
permafrost from reaching that tempera- 
ture. 

(1) Numerical analysis of a number 
of pile installations indicates that pile 
spacing should be at least five diam- 
eters of the drill hole size. A plot, similar 
to that shown in figure 4-4, may be 
prepared to relate pile spacing and 



permafrost temperature rise for the 
volumetric latent heat of the slurry 
backfill introduced into the drill hole. 
A family of curves may be developed to 
account for variation of slur ry volumet- 
ric latent heat. 

(2) In this example the slurry back- 
fill was placed at a temperature slightly 
above freezing (33.5°F) and, theoreti- 
cally, the sensible heat of the slurry 
should be considered. The volumetric 
capacity of the unfrozen slurry was 
[72(0.17 + 1.0 x 0.45) =] 44.6 Btu/ft 3 °F, 
and with a temperature difference of 
(33.5 - 32 =) 1.5°F, this represents a 
sensible heat of (1.5 x 44.6 =) 67 Btu/ft 3 . 
A comparison of this quantity with the 
volumetric latent heat of the slurry 
(4670 Btu/ft 3 ) shows that its heat may 
be considered negligible, as long as it 
near the freezing point. 

d. Artificial freeze-back time. If per- 
mafrost is temperatures are marginal, 
it may be necessary to refrigerate the 
pile to accelerate slurry freeze-back 
time and to have refrigeration available 
if permafrost temperatures rise after 
construction. The following example 
shows calculations required to deter- 
mine the amount of heat to be extracted 
from the ground. The average volume 
of slurry backfill for a group of piles is 
31 cubic feet each. The slurry is placed 
at a n average temperature of 48°F and 
must be frozen to 23°F. A silt-water 
slurry of 80 lb/ft 3 dry weight and 40 
percent water content is used as back- 
fill material, and an available refrigera- 
tion unit is capable of removing 225,000 
Btu/hr. Calculate the length of time 
required to freeze back a cluster of 20 
piles. 

—Volumetric latent heat of back- 
fill: 

L = (144 x 80 x 0.40) = 4600 Btu/ft 3 . (eq 4-16) 

—Volumetric heat capacity of 
frozen backfill: 

C f = 80 [0.17 + (0.5X0.4)] = 29.6 Btu/ft 3 »F. 

(eq 4-17) 

—Volumetric heat capacity of un- 
frozen backfill: 

C u =■ 80 (0.17 + (1.0X0.4)] = 45.6 Btu/ft 3 °F. 

(eq 4-18) 

— Heatrequired to lower the slurry 
temperature to the freezing point: 

45.6 X 31 (48 - 32) = 22,618 Btu/pile. (eq 4-19) 

—Heat required to freeze slurry: 

31 X 4600 » 142,600 Btu/pile. (eq 4-20) 
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— Heat required to lower the 
slurry temperature from the freezing 
point to 23°P; 

39.6 X 31 (32 ■ 23) = 8258 Btu/pile. (eq 4-21) 

—Total heat to be removed from 
the slurry: 

20 (22,600 + 142,700 + 8200) = 3,470,000 Btu. 

(eq 4-22) 

—Time required for artificial 
freeze-back, excluding allowances for 
system losses: 

3,470,000/225,000 - 15.5 hours. (eq 4-23 

e. Heat transfer by thermal piles. 
Artificial freeze-back may be accom 
plished also by use of two types of self 
refrigerated heat exchangers: a single 
phase liquid -convection heat transfer 
device and a two-phase boiling-liquid 
and vapor convection heat transfer 
device. TM 5-852-4/ AFM 88-19, Chapter 
4 presents heat transfer rates for the 
two-phase system. There are few heat 
transfer field data available for the 
single-phase system. 

4-3. Utility distribution systems in frozen 

ground. 

General considerations for the design 

of utility systems in cold regions are 

given in TM 5-852 -5/AFR 88-19, Volume 

5. 

a. Burying ivater pipes in frozen 
ground. Water pipes that are buried in 
frozen ground may be kept from freez- 
ing by any one of the following 
methods: 1) placing the water line in an 
insulated utilidor, which is a contin- 
uous closed conduit with all lines, such 
as water, sewage and steamlines, in- 
stalled away from direct contact with 
frozen ground, 2) providing a suf- 
ficient flow velocity such that the water 
temperature at the terminus of the pipe- 
line does not reach the freezing point 
or 3) heating the water at the intake or 
at intermediate stations along the line. 
A layer of insulation around a pipeline 
will retard, but not prevent, freezing of 
standing water in a pipe. The thermal 
analysis of a pipeline buried in frozen 
ground is complicated by the changing 
thermal properties, ice content, sea- 
sonal and diurnal changes of tempera- 
ture and the intermittent water flow. 
Some calculation techniques applic- 
able to the problem of buried utilities 
are presented below for standing and 
for flowing water. Additional tech- 



niques are presented in TM 5-852- 
5/AFR 88-19, Volume 5. 

b. Freezing of standing water in a 
buried pipe. Problems with freezeup 
of stationary water must take into 
account the initial time required to 
lower the water temperature to the 
freezing point and the amount of time 
required to form an annulus of ice in 
the pipe. In most instances the danger 
point is reached when water begins to 
freeze. 

(1) The time required to lower the 
temperature of nonflowing water in 
an insulated pipe to the freezing point 
is given by the equation 



31.2 



K ; 



0-2 In 



. -T = 



)ln 



32 -TV 



(eq 4-24) 



where 

t 



p 
**i 

T„ 



= time (hr) 

= thermal conductivity of in- 
sulation (Btu/ft hr °F) 
= radius of pipe (ft) 
= radius to outer edge of insula- 
tion (ft) 
= initial water temperature (°F) 
T s = temperature of surrounding 
soil (°F). 
For example, a 12-inch diameter iron 
pipe containing water at42°F is buried 
in 28°F soil. Determine the time re- 
quired to lower the water temperature 
to 32°F if the pipe is insulated with 3 
inches of cellular glass (K, = 0.033 
Btu/ft hr°F). 



31.2 
0.033 



(0.50' 



0.75 42-28 

In m ) In 



0.50 



32-28 



(2) 



120 hours (5 days) 
eq 4-25 

Once the water temperature 



has been lowered to the freezing point, 
ice begins to form in an annular ring 
inside the pipe. The following assump- 
tions are made to solve this problem: 

—The water is initially at 32°F. 

— The heat released by the freezing 
of water does not affect the sur- 
rounding ground temperatures. 

— The volumetric heat capacity of 
the ice may be ignored. 

—The thermal resistance of the 
pipe wall is negligible. 
The solution predicts the time required 
to form an annulus of ice around the 
inner wall of the pipe. Knowledge of 
pipe radius, insulation thickness and 
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and thermal properties, thermal con- 
ductivity of ice, latent heat of fusion of 
water, and surrounding ground tem- 
peratures are necessary to solve this 
problem. The temperature of ground 
surrounding the pipe and the time 
during which the ground remains 
below freezing is difficult to estimate. 
The relationship between time and the 
radius of ice formed inside an insulated 
pipe is given by the expression 

i-i r 2 v 2 

in — -+ 1/2X1- ~o-)-{ — -)ln -H. 
P r p P 

(cci 4-20) 



LrJ 



2KAT 



[( 



K 



where 
t 
L 



r P 
K 



= time (hr) 

= latent heat of water (9000 

Btu/ft 3 ) 
= radius of pipe (ft) 

= thermal conductivity of ice 
(1.33 Btu/ft hr °P) 

AT = temperature difference be- 
tween water and surround- 
ing soil (°F, assume water tem- 
perature is 32°F) 

K, = thermal conductivity of in- 
sulation (Btu/ft hr °F) 

r, = radius to outer edge of insula- 
tion (ft) 

r = inner radius of ice annulus 
(ft). 
If the pipe is not protected by insula- 
tion, the equation is 



Lrf 



2KAT 



[1/2(1- 



)-( 



) 2 ln 



]. (eq4-27) 



This expression for an insulated pipe 
may be simplified by rearrangement 
and substitution of numerical values 
for the latent heat and the thermal con- 
ductivity of ice. This yields 



' 1690 — ly 



where 



[i- 



<l-ln 



)]• 



(eq 4-28) 



(eq 4-20) 



r P P 

The relationship between y and r/r p is 
given in figure 4-5. 

Following is an example. A 12-inch iron 
pipe, insulated with 3 inches of cellular 
glass (Kj ■ 0.033 Btu/ft hr °F), is placed 
in 28°F soil. Calculate the time required 
to reduce the bore of the pipe to 6 



inches and the time required to com- 
pletely freeze the water. Assume the 
rate of flow does not influence freezing. 
The time required to reduce the bore 
to 6 inches will be 



OOOO (0.5) 3 
2<1.33)(32-38) 



1.33 



0.75 



[( In +0.5) 

0.033 0.50 



0J35 2 0.25 2 0.50 

(1 -d^ 2) - ( oTo 1 ln 6^ ] 

» 2640 hours. 



(eq 4-30) 



The time required to completely freeze 
the water in the pipe will be 



t - 



OOOO (0.5) 2 

2(1.33X32-38) 




1.33 



0.75 



[( In + 0.5) 

0.033 0.50 



0.50 

(1 -o^ 2) - ( oTo )2ln -o-l 

■ 3550 hours. (eq4-31) 

The calculation is simplified since the 
term "r" for the inner radius of the pipe 
goes to zero. For an uninsulated pipe, 
the calculations assume the simplified 
form of 



t - 



(1690)(0.5) 2 
(32-28) 



(eq 4-32) 



where r/r p - and y = 1.0 (fig. 4-5). 
Therefore, 

t = 106 hours. 
This example illustrates the effective- 
ness of insulation in retarding the 
freezeup of water in pipes, but as stated 
above, the assumptions used to develop 
these equations are conservative and 
the actual length of freezing time would 
be greater. 

c. Thawing of frozen soil around a 
suddenly warmed pipe. In the pre- 
ceding example it was assumed that 
water initially at 32°F was placed in 
frozen ground and the time relation- 
ship for freezing of the water in the 
pipe was determined. If the water was 
maintained above freezing, the frozen 
soil surrounding the ppe would thaw. 
To formulate a mathematical expres- 
sion relating the time with the radius 
of thaw, it is assumed that: 1) the volu- 
metric heat capacity of the soil is neg- 
ligible, 2) both the surrounding soil 
and pipe are initially at32°F, and 3) the 
pipe temperature is suddenly raised to 
a temperature above 32°F. The formula 
for an insulated pipe is 
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Figure 4-5. Freezeup of stationary water in an uninsulated pipe. 



Lr* 



[( 



K. 



K« 



[1- 



3T > +ln 



In 



] 



•0.5) 



(eq 4-33) 



where 

t ~ time (hr) 

L - latent heat of soil (Btu/f t 3 ) 
r = radius to outer edge of thaw- 
ing soil (ft) 



K u = thermal conductivity of un- 
frozen soil (Btu/ft hr °F) 

AT = temperature difference be- 
tween pipe and surrounding 
soil (°F, assume soil at 32°F) 

K, = thermal conductivity of in- 
sulation (Btu/ft hr °F) 

r i » radius to outer edge of insula- 
tion (ft) 

r p = radius of pipe (ft). 
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If the pipe is not protected, by insula- 
tion, the expression is 



Lr< 






[-0.5 (1 ■ 



) + ln 



]. (eq4-34) 



For example, a 7- by 7-foot concrete 
utilidor has 9-inch concrete walls with 
an outer covering of 6 inches of insula- 
tion (K, - 0.033 Btu/fthr°F). The frozen 
soil around the utilidor is a sandy 
gravel with a dry density of 115 lb/ft 3 
and a water content of 7.8 percent at a 
temperature of 32°F. Neglect the ther- 
mal resistance of the concrete, and 
determine the time required to thaw 1 
foot of soil when the temperature of 
the utilidor walls is suddenly raised to 
50°F. 

K u = 1.2 Btu/ft hr °F 

L = 1290 Btu/ft 3 

K t = 0.033 Btu/ft hr °F. 
For calculation, square sections may 
be treated as cylinders of the same 
perimeter. 



Symbol 



Dimension-square 
(ft) 



Equivalent radius 
(ft) 



r 9.5 + 24/12=11.5 

r, 8.5+ 12/12= 9.5 

r„ 7.0 + 18/12 = 8.5 



7.33 
6.04 
5.42 



t» 



<1290)(7.33) 3 
3U.2)<50-32) 



(1- 



6.04^ 



7.33-* 



1.2 



0.033 



7.33 



In 



6.04 



5.42 



■0.5) 



) + ln ] 

' 6.04 



= 2060 hours. 



(eq 4-35) 



d. Practical considerations. 

(1) The above-mentioned formulas 
indicate the relationship between time, 
radius of freeze or thaw, and tempera- 
ture difference between the water in 
the pipe and the ground. In sufficient 
time, standing water in the pipe will 
freeze or frozen ground will thaw, 
depending on temperature differen- 
tials. For practical problems the as- 
sumed constant temperature differen- 
tials will not exist for a long time but 
will vary with season and even with 
the hour at shallow depths. The freez- 
ing and thawing index concept con- 
siders the intensity of temperature dif- 
ferential from freezing (32°F) and the 
duration of this differential. In the 



preceding equations, the time t can be 
multiplied by the temperature differen- 
tial AT to give either a freezing or 
thawing index, and the radius of ice 
formation or thawed ground radius 
can be determined by trial-and-error. 
It was shown that a temperature dif- 
ferential of (32-28) = 4°F lasting for 3150 
hours would result in complete freeze- 
up of the water in the pipe; this is equiva- 
lent to a freezing index of [(4 X 3150)/24 
=] 524 degree-days. If the freezing 
index at the depth of pipe burial were 
less than 524, the standing water in the 
pipe would not completely freeze in 
that time. Thus, the freezing index at a 
particular depth can be used to fore- 
cast freezeup of stationary water in 
pipes located in the annual frost zone. 
(2) Even insulated water lines 
located in frozen ground usually re- 
quire an inlet water temperature sig- 
nificantly above freezing. Whether the 
water lines are insulated or not, this 
may thaw some of the surrounding 
frozen ground. This thawed annulus 
will retard water freezeup in the pipe 
if and when flow conditions change. 
The situation is practice is generally 
complicated by the intermittentcharac- 
ter of water demand. In some northern 
communities the problems of irregular 
water demand are solved by con- 
structing the water lines in a contin- 
uous loop with provisions for periodic 
flow reversals. Water temperatures 
should be closely monitored and water 
usage patterns considered in esti- 
mating water freezeup. 

e. Freezing and flowing water in 
buried pipes. Problems involving freez- 
ing of flowing water in buried pipes 
require knowledge of the distance the 
water will flow before the temperature 
of the water lowers to the freezing 
point. By providing enough above- 
freezing water, the loss of heat to the 
surrounding frozen soil can be 
balanced to provide an outlet tempera- 
ture slightly above freezing. The prob- 
lems of freezing of flowing water in 
insulated and bare pipe are illustrated 
below. 

(1) Insulated pipe. It is assumed 
that 1) the temperature of the frozen 
ground surrounding the pipe is con- 
stant for the period of flow over the 
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entire length of pipe, 2) the effect of 
friction heat developed by water flow 
is negligible, 3) the thermal resistance 
and heat capacity of the pipe wall are 
negligible, and 4) the temperature dis- 
tribution of the water in the pipe is 
uniform at each cross-sectional area. 
The velocity required to prevent freeze- 
up of flowing water in a pipe is given 
by 



sK, 



2 Fi 
112,000 r~ {In )(ln 

P r P 



T 






) (eq4-36) 



where 

V = velocity of flow (ft/s) 
s = length of pipeline (ft) 
K } = thermal conductivity of in- 
sulation (Btu/ft hr °F) 
r p = radius of pipe (ft) 
r t = radius to outer edge of insula- 
tion (ft) 
Tj = inlet water temperature (°F) 
T 2 = outlet water temperature (°F) 
T s = temperature of surrounding 
frozen soil ( C F). 
For example, an 11,000-ft long, 6-inch- 
diameter pipe is buried in 10°F soil. 
The pipe is covered with a 2-inch layer 
of insulation (Kj = 0.03 Btu/ft hr°F) and 
the inlet water temperature is 39°F. 
Calculate the velocity of flow required 
to keep the water from freezing. 

(11,000)(0.03) 



v 



0.417 39 - 10 

112,000 <0.25) 2 (ln )(ln ) 

~~" 32 -lO 



0.25 
= 0.33 ft/s (20 ft/min). (eq 4-37) 

To provide for temporary reductions 
in flow and in recognition of the un- 
certainties concerning the manner of 
ice formation within the pipe, it is recom- 
mended that the velocity of flow be 
doubled in design. 

(2) Uninsulated pipe. 
(a) When flowing water is first 
introduced into a bare pipe buried in 
frozen ground, the heat loss from water 
is greater than it is after the system 
has been in operation for a period of 
time. The initial heat loss is greater 
because the pipe wall and the soil im- 
mediately adjacent to the pipe are 
colder than they are after water has 
flowed over a time. Soil temperatures 
surrounding the pipe increase and 



eventually become reasonably stable 
with time. An expression relating these 
variables is 



! exp 



2iv 



h 

v 



5.6X10 4 



(eq 4-38) 



where 

Tj = inlet water temperature (°F) 
T s = frozen soil temperature (°F) 
T 3 = ou tlet water temperature (°F) 
s = length of pipeline (ft) 
2r p = diameter of pipe (ft) 
h =heat transfer coefficient 

(Btu/ft 3 hr °F) 
V = velocity of flow (ft/s). 

A nomogram of this equation is shown 

in figure 4-6. 

(b) Limited field experiments in 
clay and sandy clay soils suggest 
values of h for metal pipelines subject 
to normal use (conditions or intermit- 
tent flow) of 6.0 for the initial period of 
operation an 2.0 thereafter. These 
values are not applicable for pipes 
smaller than 4-inches in diameter. The 
h value is dependent upon the thermal 
properties of the surrounding soil, the 
diameter of the pipe, the type of pipe 
material and the temperature gradient 
in the ground around the pipe's radius. 
The value for h given above provides a 
reasonable basis for design of pipelines 
in which the total quantity of water 
consumed per day is at least eight times 
the volume of pipes in the entire system. 
The time of operation required for the 
temperature distribution in the water 
to stabilize is approximately. 



■ 0.005 ■ 



v 



(eq4-39) 



where 

S 

V 



= time (hr) 

= length of pipeline (ft) 

= velocity of flow (ft/s). 



fcjFor example, water atan inlet 
temperature of 40°F flows at 2 ft/s in a 
12-inch iron pipeline, 2.2 miles long. 
The ground temperature surrounding 
the pipe is 25°F. Estimate the outlet 
water temperature during the initial 
period of flow (h = 6.0) and after sur- 
rounding temperatures have stabilized 
(h - 2.0). 
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— Initial period (see fig. 4-6): 

s 2.2 X 5280 



2r, 



- U6X10 4 (eq 4-40) 



V =■ 2 ft/s. 

h ■ 6.0 Btu/ft 2 hr °F 



thus 



L l -Tc 



1.77 ' 



T 2 T S ' 
To - 33.5°F. 



40-25 
To -25 



(eq 4-41) 



—Stabilized period: 

s 
2i\ 



1.16X10 4 

2 ft/s 
h= 2.0 Btu/ft 2 hr °F 



P 

V » 2 ft/s 



thus 



T, -Tc 



' 1.22 



.40 - 25 
To -25 



(eq 4-42) 



Tr, - 37.3°F. 



Figure 4-6. Temperature drop of flowing water In a pipeline. 
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To operate on the safe side, the outlet 
temperature T 2 should remain at or 
above 35°F. The calculated initial water 
temperature of 33.5°F would be con- 
sidered unsafe and either the flow velo- 
city should be increased to approxi- 
mately 3 ft/s or the inlet water tempera- 
ture should be raised to about 43°F. 
These^ precautions would be needed 
only for an initial period, i.e., 



11,600 



t„ = 0.005 — = 0.005 
° V 



= approximately 29 hours. <eq 4-43) 

/. Design considerations. The cal- 
culation techniques presented above 
indicate the principal factors to be con- 
sidered in design of water distribution 
lines placed in frozen ground. The use 
of these techniques together with recog- 
nition of the complexities of actual in- 
pipe ice formation and sound engi- 
neering judgment provides a basis for 
design of pipelines in areas of seasonal 
frost and permafrost Changing the sur- 
face cover over an installed pipeline 
will affect the distribution of tempera- 
tures with depth and may result in 
depressing the temperatures adjacent 
to the pipe. This is particularly true if a 
natural vegetative cover is stripped and 
replaced by asnow-free pavement. The 
influence of new construction above 
an existing pipeline may require a 
change in operating procedures for 
the system, such as an increase in the 
velocity or flow or additional heating 
at the inlet. 

4-4. Discussion of multidimensional heat 
flow. 

a. The relatively simple analytical 
techniques discussed in this manual 
are not always sufficient for consider- 
ing the concurrent thermal effects of 
multidimensional temperature change 
and soil water phase transformation. 
The one-dimensional and radial heat- 
flow computation techniques pres- 
ented in this manual were based on 
field observations and the use of reason- 
able simplifying assumptions. The tech- 
niques are intended to facilitate ana- 



lysis and to promote adequate design. 
The assumptions involved and tech- 
nique limitations have been em- 
phasized. 

b. Heat flow beneath heated struc- 
tures is multidimensional because of 
the finite boundaries of such struc- 
tures. The ground surface temperature 
adjacent to the south side of the build- 
ing is generally higher than that on 
the north side, and the ground surface 
temperature on the west side is gen- 
erally higher than on the east side; in 
addition to this influence, the three- 
dimensional temperature distribution 
beneath the building will be affected 
by the plan dimensions of the floor. 
Three-dimensional solutions are avail- 
able to the problem of heat flow in 
homogeneous materials beneath the 
surface of a heated finite area sur- 
rounded by an infinite area subject to 
adissimilar surface temperature condi- 
tion; however, the solutions consider 
only the effects of temperature change 
and not the effects of phase transforma- 
tions. Such solutions tend to be rather 
complex and unwieldy, and their 
neglect of latent heat generally results 
in an over-estimation of the depth of 
freeze or thaw. The magnitude of this 
over-estimation is dependent on the 
quantity of moisture in the frozen or 
thawed soil. 

c. The example given in paragraph 
3-8a for calculating the depth of thaw 
beneath a heated slab-on-grade build- 
ing considered only one-dimensional 
vertical heat f low and excluded lateral 
heat flow from the soil beneath the 
building to the surrounding soil mass 
in the winter. The amount of lateral 
heat flow would depend on the building 
dimensions and the wintertime soil tem- 
perature gradient. Slab-on-grade, 
heated structures usually prevent 
frost penetration under the center of 
the building and result in a thaw bulb 
in the foundation soil that may cause 
permafrost degradation with time. This 
type of construction is discussed in 
TM 5-852-4/ AFM 88-19, Chapter 4. 
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APPENDIX A 



Symbol 

a 
a' 
A 

A d 
c 

c p 
C 

c f 



d 

e 
erf 



exp(x) 

f 

F 

f 

h 

h f 

h v 
H 

I 

If 
k 

K 

K f 

K u 

L 
m 
In 
I 



Description 

Thermal diffusivity 

Shortest dimension 

Sinusoidal temperature amplitude 

Duct cross-sectional area 

Specific heat 

Specific heat of air at constant pressure 

Volumetric heat capacity 

Volumetric heat capacity in frozen 

condition 

Volumetric heat capacity in 
unfrozen condition 
Thickness of soil layer 
Equivalent duct diameter 
Roughness factor 
Error function; 
erf z ■ (2/iH) o^e-Ai/i, 
where erf °° = 1 and erf (-z) 
= -erf z 

e x 

Frequency of sine wave 

Air freezing index 

Friction factor 

Heat transfer coefficient 

Surface conductance for combined 

radiation and convection 

Draft head 

Friction head 

Velocity head 

Stack height 

Air thawing index 

Floor thawing index 

Coefficient of thermal conductivity 

Thermal conductivity 

Thermal conductivity in frozen 

condition 

Thermal conductivity in unfrozen 

condition 

Volumetric latent heat of fusion 

Duct spacing 

Natural logarithm 

Duct length 



Units 



ft 2 /day or ft 2 /hr 


ft 


op 


ft 2 


Btu/lb °F 


Btu/lb °F 


Btu/ft 3 °F 


Btu/ft 30 F 


Btu/ft 3 °F 


ft 


ft 


ft 


dimensionless 



dimensionless 
cycles/day 
degree-days 
dimensionless 
Btu/ft 2 hr °F 

Btu/f t 2 hr °F 
inches of water 
inches of water 
inches of water 
ft 

degree-days 
degree-days 
Btu/fta hr °F per in. 
Btu/ft hr °F 

Btu/ft hr °F 

Btu/ft hr °F 

Btu/ft 3 

ft 

dimensionless 

ft 
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Symbol Description 

ib 

to 

is 

MAT 



n 

nP 

nl 

N 

P 

Q 

r 

R 

S 

s 

t 

T 

T 



Units 



R 



R 



T 

T s 

T s 

T„ 

V 

v 

V s 
W 

X 

a 
Ya 

£ 

v 

IT 
P 




Allowance for bends, etc., in duct 


ft 


Equivalent duct length 


ft 


Length of straight duct 


ft 


Mean Annual temperature 


op 


surface index 




"n"-factor » . . . 

air index 


dimensionless 


Surface freezing index 


degree-days 


Surface thawing index 


degree-days 


Reynolds number 


dimensionless 


Period of sine wave 


365 days 


Latent heat per linear foot of slurry backfill 


Btu/ft 


Radius 


ft 


Thermal resistance 


ft s hr °P/Btu 


Pile spacing 


ft 


Length of pipeline 


ft 


Time 


hours or days 


Temperature 


°p 


Initial temperature 


°F 


Temperature rise in duct air 


°F 


Surface temperature 


op 


Soil temperature 


op 


Water temperature 


op 


Velocity of flow 


ft/s 


Initial temperature of the soil with respect 




to32°F 


op 


Average surface temperature with respect 




to32°P 


op 


Water content 


percent dry weight 


Depth of freeze or thaw 


ft 


Thermal ratio ■ v /v s 


dimensionless 


Dry unit weight 


lb/ft 3 


Efficiency 


percent 


Lambda coefficient 


dimensionless 


Fusion parameter 


dimensionless 


Kinematic viscosity 


fts/hr 


Pi 


3.14 


Density of air 


lb/ft 3 


Heat flow to duct 


Btu/ft 2 per hr 
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APPENDIX B 
THERMAL MODELS FOR COMPUTING FREEZE AND THAW DEPTHS' 



B- 1 . Analytical Solutions. 

a. The Oxford University Press 
publication, Conduction of Heat in 
Solids, 2 presents solutions to many 
one-, two- and three-dimensional heat 
flow problems. Homogeneous isotropic 
materials are used in most solutions, 
but some solutions are presented for 
layered systems. U.S.G.S. Bulletin 1083- 
A 3 developed a one-dimensional tech- 
nique for predicting the damping of a 
periodic surface temperature at dif- 
ferent depths in a two- and three- 
layered soil system. Unidirectional 
heat flux was considered. U.S.G.S. Bul- 
letin 1052-B 4 developed a method for 
estimating the three-dimensional ther- 
mal regime in a homogeneous isotropic 
soil beneath a heated structure. None 
of these techniques considers phase 
change of the soil moisture. Neglect of 
the effects of latent heat of fusion of 
the soil moisture normally does not 
cause substantial error in prediction 
of frost depths in soil of low water 
content. Differences between actual 
and computed thaw depths increase 
rapidly with increasing water content 
because of the increased volumetric 
heat capacity and greater latent heat 
of the wetter soil . 

b. Several empr ileal and semi-empir- 
ical equations have been developed that 
consider latent heat of fusion of the 
soil, the Stefan equation was originally 
developed to calculate the thickness of 
ice on a calm body of water (isothermal 
at the freezing temperature, 32°F) 
expressed as; 

X, =V 48K i F/L i (eqB-1) 



The documents mentioned in this appendix 
are sources for additional information and 
are found in the bibliography. 

2 Carslaw and Jaeger, 1959. 

3 Lachenbruch, 1959. 

4 Lachenbruch, 1957. 



where 

Xj = ice thickness (ft) 
K t » thermal conductivity of ice 
(Btu/ft hr °P) 

F = freezing index (degree-days) 
Lj = volumetric latent heat of 
fusion of ice (Btu/ft 3 ). 
The Stefan equation has been modified 
by many individuals and agencies, and 
many similar equations have been 
developed. Some of the equations use 
functions or initial conditions slightly 
different from those used in the 
original Stefan model. The most widely 
used equation to estimate seasonal 
freeze and thaw depths is the modified 
Berggren equation. Application of this 
equation, given below, has been very 
widespread in North America. 
USACRREL Special Report 122 5 
developed a computer program for cal- 
culating freeze and thaw depths in 
layered systems u sing this equation: 
X=y48KnF/L or X = /\y48KnI/L 

(eq B-2) 
where 

X - depth of freeze or thaw (ft) 

K - thermal conductivity of soil 

(Btu/ft hr °F) 

L » volumetric latent heat of fusion 
(Btu/ft 3 ) 

n = conversion factor from air in- 
dex to surface index (dimen- 
sionless) 

F =air freezing index (degree- 
days) 

I =air thawing index (degree- 
days) 

^ = coefficient that considers the 
effect of temperature changes 
within the soil mass. It is a 
function of the freezing (or 
thawing) index, the mean 
annual temperature and the 
thermal properties of the soils. 



Aitken and Berg, 1968. 
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c. An equation very similar to the 
Stefan equation is used in the USSR to 
calculate the "standard" depth of freez- 
ing for foundation designs. Many other 
closed-form analytical techniques are 
also used in the USSR. 

B-2. Graphical and analog methods. 

a. Graphical methods have been 
used to estimate depths of freeze and 
thaw. The flow net technique can be 
used to estimate steady-state tempera- 
ture conditions. National Research 
Council of Canada, Technical Paper 
No. 163, 6 presents a graphical means to 
determine temperature in the ground 
under and around natural and engi- 
neering structures lying directly on 
the ground surface. 

fc>. Analog techniques are also used 
to estimate freeze and thaw depths. 
Table B-l shows thermal, fluid and elec- 
trical analogies. Electric analog com- 
puters are available, cost relatively 
little and are reasonably simple to use. 
The primary disadvantage of these 
machines are that reprogramming is 
normally necessary for each problem 
and complex geometries are difficult 
to simulate adequately. Hydraulic 
analogs are also available. The primary 
disadvantages of these computers are 
their complex tubing systems, space 
requirements, and the necessity to 
thoroughly clean and reconstruct them 
for each problem. At any instant of 
time, however, hydraulic analogs 
graphically show the temperature 
distribution. 



B-3. Numerical techniques. 

a. Because of the widespread avail- 
ability of electronic digital computers, 
their application to numerical solu- 
tions of the continuity equation is com- 
monplace. Numerical procedures are 
approximations to the partial differen- 
tial equation; however, they are much 
more accurate and versatile in solving 
complex transient heat flow problems 
than are the analytical techniques. 
Many computer programs allow flex- 
ible definition of boundary and initial 
conditions for both one- and two- 
dimensional problems. General back- 
ground information on the finite dif- 
ference methods available to solve heat 
flow problems are discussed in Inter- 
national Textbook Company's Heat 
Transfer Calculations by Finite 
Differences. 7 Since rectangular ele- 
ments are normally used, complex 
geometries are difficult to simulate ac- 
curately unless small element sizes are 
used. 

b. Use of the finite element technique 
is also widespread (see The Finite 
Element Method in Engineering 
Science 8 or The Finite Element Method 
in Structural and Continuum 
Mechanics,® both from McGraw-Hill). 
Elements of various shapes can be used 
with this technique; but the triangular 
shape is commonly used for two- 
dimensional problems. Complexbound- 
ary geometries can be more closely 



6 Brown, 1963. 



Dusinberre, 1961. 
Zienkiewicz, 1967. 
Zlenkiewlcz, 1971. 



Table B-1. Thermal, fluid and electric analogs (U.S. Army Corps of Engineers). 













Medium 






Item 


Thermal 




Fluid 




Electric 




A - Variables 


(D 
(2) 
(3) 


Heat 
Heat flux 
Temperature 


fJ 

q 

T 


Volume 

Flow 

Head 


S 
Q 

h" 


Charge density 
Current density 
Voltage 


J> 

T 

e 



B - Principles: 






Continuity 


(D 


9 " -» - 


i V ■ q 






9t 


Conductivity 


(2) 


q = -kVT 


Capacitance 


(3) 


d/u = CdT 



3S 



+ V ■ Q = 



+ V • j = 



dt 

Q = -kVH 

dS = AdH 



at 

1 = -aVe 
pdV = Cde 
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simulated using finite element pro- 
cedures. For multidimensional heat 
flow problems, the finite element pro- 
cedure is frequently more efficient, i.e., 
it requires less computer time than the 
finite difference technique. 

c. Many flexible computer programs 
exist that simulate heatconductionand 
phase change in soils. Each has its 
own particular datarequirements,com- 
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putational capabilities, and acquisition 
costs and restraints. USACRREL has 
completed and documented a model 
that may be useful in solving many 
analytical problems related toconstruc- 
tion in the Arctic and Subarctic; other 
models are under development. Con- 
tact HQ (DAEN-ECE-G) or HQ AFESC 
for assistance in selecting an appro- 
priate model. 
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APPENDIX C 
REFERENCES 



Government Publications 

Departments of the Army and the Air Force. 

TM 5-818-2/ AFM 88-6, Chapter 4 Pavement Design for Seasonal Frost 

Conditions 
TM 5-852 -1/AFR 88-19, Volume 1 Arctic and Subarctic Construction, 

General Provisions 
TM 5-852-4/ AFM 88-19, Chapter 4 Arctic and Subarctic Construction, 

Building Foundations 
TM 5-852-5/ AFR 88-19, Volume 5 Arctic and Subarctic Construction, 

Utilities 

Nongovernment Publications 

American Society of Heating, Refrigeration, and Air-Conditioning Engineers 

(ASHRAE) 

345 East 47th Street, New York, 

N.Y. 10017 

ASHRAE Guide and Data Book (1963) 
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